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Research  in  convective  heat  transfer  on  internal  separated  flows  has  been  extensively  conducted  in  the 
past  decades.  This  review  summarizes  numerous  researches  on  two  topics.  The  first  section  focuses  on 
studying  the  fluid  flow  and  heat  transfer  behavior  of  different  types  of  single-phase  fluid  flows  over  back¬ 
ward  facing  step  (BFS)  at  different  orientations.  The  second  section  concentrates  on  everything  related 
to  nanofluids;  its  preparation,  properties,  behavior,  applications,  and  many  others.  The  purpose  of  this 
article  is  to  get  a  clear  view  and  detailed  summary  of  the  influence  of  several  parameters  such  as  the  geo¬ 
metrical  specifications,  boundary  conditions,  type  of  fluids,  and  inclination  angle  on  the  hydrodynamic 
and  thermal  characteristics  using  (BFS).  The  reattachment  length  and  maximum  Nusselt  number  are  the 
main  target  of  such  research  where  correlation  equations  were  developed  and  reported  in  experimental 
and  numerical  studies.  The  heat  transfer  enhancement  of  nanofluids  along  with  the  nanofluids  prepara¬ 
tion  technique,  types  and  shapes  of  nanoparticles,  base  fluids  and  additives,  transport  mechanisms,  and 
stability  of  the  suspension  are  also  discussed. 
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1.  Introduction 

The  phenomena  of  flow  separation  and  subsequent  reattach¬ 
ment  which  occur  due  to  a  sudden  expansion  in  the  flow  passages 
such  as  backward  facing  steps  have  been  recognized  as  impor¬ 
tant  industrial  situations.  This  complex  flow  structure  present  in 
heating  or  cooling  applications  such  as  cooling  electronic  equip¬ 
ments,  cooling  turbines  blades,  combustion  chambers,  chemical 
processes,  cooling  of  nuclear  reactors,  wide  angle  diffusers,  high 
performance  heat  exchangers,  energy  systems  equipment,  and  flow 
in  valves.  In  many  instances  separation  of  flow  is  undesirable  and 
leads  to  unwanted  pressure  drops  and  energy  losses  which  require 
additional  fan  or  pumping  power  to  overcome  them.  However,  in 
other  circumstances  flow  separation  may  be  encouraged,  such  as 
in  burner  flame  stabilization  use  for  turbulence  promotion  leading 
to  enhanced  mixing  or  heat  and  mass  transfer  rates. 

The  problem  of  flow  over  backward  facing  step  has  been 
extensively  investigated,  both  experimentally  and  numerically. 
However,  case  studies  such  as  injection  flow,  ribs  geometry,  double 
step  ducts,  and  onset  inlet  flow,  are  not  considered  in  the  present 
review.  The  current  review  will  only  consider  the  studies  of  the 
case  of  steady-state  convective  flow  and  heat  transfer  over  single 
backward  facing  steps. 

A  large  number  of  experimental  and  numerical  studies  focus 
on  flow  and  heat  transfer  behavior  of  convective  flow  over  back¬ 
ward  facing  step  geometry  have  been  reported.  These  vast  results 
and  information  have  dealt  with  different  conditions,  parameters, 
geometry  dimensions,  and  instrumentation,  which  indeed  to  unde¬ 
fined  solid  base  for  comparison  purposes  to  indicate  more  accurate 
methodology  for  solving  the  case  studies. 

The  dispersion  of  these  results  were  attracted  the  attention  of 
many  researchers  to  unify  the  information  to  a  general  criteria. 
In  1992  and  1993,  the  Aerospace  Heat  Transfer  Committee  (K- 
12)  of  the  Heat  Transfer  Division  of  the  ASME  held  a  technical 
session  for  benchmark  heat  transfer  problems  during  the  Winter 
Annual  Meeting  to  solve  numerically  laminar  mixed  convection 
flow  over  two-dimensional  horizontal  and  vertical  backward  fac¬ 
ing  step,  respectively  [1-3].  Moreover,  due  to  the  importance  of  the 
separation  and  reattachment  phenomenon  Abu-Mulaweh  in  2003 
[4]  reviewed  the  results  of  flow  and  heat  transfer  of  single-phase 
laminar  mixed  convection  flow  over  different  orientations  of  both 
backward  and  forward  facing  steps  for  several  previous  work.  Most 
of  the  studies  in  this  literature  declared  that  only  one  large  primary 
recirculation  region  develops  downstream  of  the  backward-facing 
step,  while  in  some  cases  a  secondary  recirculation  may  develop 
between  the  downstream  stepped  wall  and  the  step  or  at  the  top 
wall  of  the  duct  due  to  increase  in  step  height  and  inlet  velocity. 

The  purpose  of  this  paper  is  threefold.  The  primary  purpose  is  to 
review  the  available  results  of  the  flow  and  heat  transfer  presented 
in  previously  published  data  for  single  backward  facing  step.  The 
results  of  massive  number  of  different  parameters  affect  the  fluid 
flow  and  heat  transfer  characteristics  are  summarized  and  pre¬ 
sented  with  a  large  number  of  correlation  equations  to  predict  the 
peaks  of  the  flow  and  heat  transfer  in  the  recirculation  regions.  The 
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AR  aspect  ratio  (W/s) 
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Subscripts 
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second  purpose  is  to  understand  the  characteristics  and  functions 
of  nanofluids,  to  expect  their  effects  and  heat  transfer  enhance¬ 
ment  in  such  geometries.  The  third  purpose  is  to  open  a  gate  for 
new  researches  and  propose  suggestions  that  could  lead  to  improve 
the  ability  to  predict  separation  and  reattachment  phenomenon 


H.A.  Mohammed  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  IS  (20 1 !  )  2921-2939 


2923 


using  nanofluids.  The  review  starts  with  an  extensive  review  on 
the  flow  and  heat  transfer  in  backward  facing  step  channels.  After 
that,  a  comprehensive  review  for  nanofluids  and  its  characteristics 
is  described.  At  the  end,  the  review  concentrates  on  the  flow  and 
heat  transfer  over  a  backward  facing  step  using  nanofluids. 

2.  Flow  geometry 


flow  over  a  backward-facing  step  was  done  in  the  late  1950s.  The 
advent  of  new  instrumentation  and  the  improvement  of  the  numer¬ 
ical  codes  increase  the  number  of  new  researches  in  such  problem 
and  facilitate  the  complex  study  of  three-dimensional  flow  in  the 
separation  and  reattachment  zone.  The  horizontal,  inclined,  ver¬ 
tical  cases  were  investigated  for  different  geometrical,  boundary 
conditions  and  fluids  properties. 


The  general  feature  of  backward  facing  step  flow  has  the  sim¬ 
plest  reattachment  flow,  but  the  flow  field  is  still  very  complex. 
The  scope  of  this  review  is  to  summarize  the  previous  efforts 
that  studied  the  single-phase  convective  flow  over  two  and  three 
dimensional  horizontal  single  backward  facing  step  expansion  in 
horizontal,  inclined,  and  vertical  orientations.  A  schematic  diagram 
for  backward  facing  step  at  either  uniform  wall  temperature  (UWT) 
or  uniform  heat  flux  (UHF)  is  shown  in  Fig.  1. 

The  advent  of  new  instrumentation  such  as  laser  anemometer, 
pulsed-wire  anemometer,  laser-Doppler  measurement,  hot  wire 
probes,  etc.,  have  made  a  renaissance  for  proliferation  of  new 
researches  in  this  area  and  facilitate  the  study  of  sophisticated 
three-dimensional  flow  structures  associated  with  separation  and 
reattachment  phenomena  in  this  geometry.  In  experimental  stud¬ 
ies,  the  fluid  velocity  and  temperature  distributions  were  measured 
using  Laser-Doppler  velocitimeter  and  cold-wire  anemometer, 
respectively,  at  any  location.  In  numerical  studies,  the  step  geom¬ 
etry  and  its  boundary  conditions  were  modeled  using  Boussinesq 
approximation.  The  steady  three-dimensional  mass  conservation, 
momentum  equations,  and  energy  equation  governing  the  fluid 
flow  and  heat  transfer  problem  are  reduced  to  the  following  forms 
[5]: 

Continuity  equation: 


d(pu)  3  (pv)  3  (pw)  = 
3x  +  dy  dz 

X-momentum  equation: 
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Energy  equation: 
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There  are  different  techniques  were  used  to  solve  the  partial  dif¬ 
ferential  equations  along  with  specific  boundary  conditions,  such 
as  finite  difference  (FD),  finite  volume  (FV),  finite  element  (FE),  or 
lattice  Boltzmann  (LB).  These  methods  are  solved  by  using  differ¬ 
ent  numerical  software  codes  such  as  FLOTRAN,  FLUENT,  TEACH, 
COUGAR,  KAMELEONII,  FIDAP,  NEKTON,  CUTEFLOWS,  QMRCGS, 
etc. 


3.  Backward  facing  step  geometry 

Convective  flow  and  heat  transfer  over  a  backward  facing  step 
has  been  extensively  investigated  in  the  last  decades,  both  exper¬ 
imentally  and  numerically.  It  has  been  reported  in  the  literature 
that  the  first  effort  for  studying  the  separation  and  reattachment 


3.1.  Horizontal  flow 

One  of  the  first  investigations  dealt  with  laminar  regime  and 
subsonic  free  stream  flowing  over  a  step  was  presented  by  Gold¬ 
stein  et  al.  [6].  They  concluded,  based  on  their  experimental 
measurements,  that  in  laminar  regime  the  reattachment  point  is 
not  a  constant  value  as  in  turbulent  regime,  but  it  depends  on  the 
boundary  displacement  thickness  and  the  hydraulic  Reynolds  num¬ 
ber.  Interest  in  studying  the  reattachment  phenomena  over  the 
backward-facing  step  continued  in  the  1970s  and  1980s.  Experi¬ 
mental  and  numerical  studies  for  2D  pure  forced  convection  flow 
where  there  is  no  heat  transfer  effect  have  been  investigated  by 
Denham  and  Patrick  [7],  and  Armaly  et  al.  [8,9], 

Denham  and  Patrick  [7]  studied  the  measurements  of  2D  lam¬ 
inar  flow  over  backward-facing  step.  The  water  flow  in  a  range 
of  Reynolds  number  (50  <  Re  <  250)  passed  through  a  single  plane 
duct  expansion  with  an  expansion  ratio  and  aspect  ratio  of  ER  =  3, 
and  AR  20,  respectively.  They  indicated  that  the  general  flow  field 
bears  similarities  to  other  recirculating  flows  like  axi-symmetric 
and  the  2D  duct  sudden  expansions  but  recirculation  zone  lengths 
and  recirculated  mass  flow  rates  were  smaller  at  a  given  value  of 
Reynolds  number.  They  noticed  that  a  small  fluctuation  occurs  at 
Re  =  229  while  the  other  lower  values  remain  steady;  this  fluctu¬ 
ation  observed  to  be  the  onset  of  transition  to  turbulence  of  the 
separated  boundary  layers.  Furthermore,  it  was  found  that  there 
is  a  tendency  towards  secondary  separation  at  the  upper  wall  as 


reattachment  occurs,  this  phenomenon  was  produced  due  to  the 
inflections  of  the  velocity  profiles  and  it  was  clearer  at  the  highest 
Reynolds  number. 

The  effect  of  Reynolds  number  for  different  flow  regimes  was 
conducted  by  Armaly  et  al.  [8],  They  reported  the  measure¬ 
ments  of  velocity  distribution  and  reattachment  length  through 
downstream  of  a  single  backward-facing  step  mounted  in  a  two- 
dimensional  channel.  Air  was  used  as  a  working  medium  with 
Reynolds  number  range  of  70  <  Re  <8000  in  a  test  section  has  a 
high  aspect  ratio  of  AR  =  36.  Numerical  predictions  were  also  pre¬ 
sented  to  be  compared  with  the  experimental  results  as  the  flow 
maintained  its  two-dimensionality  in  the  experiments.  The  com¬ 
parison  showed  a  good  agreement  with  the  experimental  findings 
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up  to  a  Reynolds  number  Re ^400.  Armaly  et  al.  [9]  showed  that 
the  vortices  start  to  develop  in  the  laminar  region  at  Re  ~  800, 
but  it  becomes  very  pronounced  at  Re  =  2300,  and  damped  out  at 
Re  =  5000.  It  was  indicated  that  in  both  the  pure  laminar  (Re  <  400) 
and  the  fully  turbulent  (Re  >  6600)  flow  regimes,  the  flow  is  indeed 
two-dimensional  over  large  parts  of  the  width  of  the  test  section. 

Most  of  the  existing  research  efforts  have  been  focused  on  two- 
dimensional  situations,  though  most  of  the  realistic  applications 
are  practically  dominated  by  three  dimensional  geometries.  In  the 
early  1990s  many  researchers  tend  to  study  the  flow  field  in  three- 
dimensional  due  to  the  presence  of  certain  effects  which  are  not 
visualized  in  two-dimensional  (at  centerline  of  the  plane)  such  as 
the  side-wall  effect.  The  three-dimensional  pure  forced  convection 
flow  have  been  extensively  investigated  by  Shih  and  Ho  [10],  Chi- 
ang  et  al.  [11],  De  Brederode  and  Bradshaw  [12],  Hertzberg  and  Ho 
[13],  Tylli  et  al.  [14],  Armaly  et  al.  [15],  and  Nie  and  Armaly  [16]. 
Shih  and  Ho  [10]  examined  the  water  flow  field  behind  a  small 
aspect  ratio  backward-facing  step  with  AR  =  3  using  Laser-Doppler 
Anemometer.  The  measurements  of  the  three  velocity  components 
inside  the  separation  region  indicated  that  the  flow  is  highly  three- 
dimensional.  The  mean  reattachment  line  is  highly  distorted  across 
the  span  with  the  shortest  reattachment  position  located  at  the  cen¬ 
treline  and  about  20%  shorter  than  the  value  closer  to  the  side  wall. 
The  maximum  spanwise  velocity  component  was  indicated  to  be 
of  the  same  order  as  the  streamwise  velocity  at  the  point  near  the 
center  of  the  step. 

Chiang  et  al.  [11]  aimed  to  clarify  and  present  a  deep  expla¬ 
nation  on  skin-friction  fields  as  a  guide  for  exploring  the  physical 
interpretation  of  the  three-dimensional  flow  separation  behind  the 
step,  reattachment  of  the  flow,  and  the  roof  recirculatory  flow  pat¬ 
tern.  The  channel  designed  to  have  an  expansion  ratio  ER=  1.9432 
where  the  step  was  0.9432.  The  parametric  study  was  carried  out 
by  varying  Reynolds  numbers  and  various  spans  with  values  not 


greater  than  800  and  10,  respectively.  It  was  shown  that  the  flow 
separation-attachment  occurs  on  the  roof  of  the  channel  but  it  con¬ 
fined  only  to  the  sidewall  region.  The  roof  recirculation  zone  size  is 
dependent  on  Reynolds  number. 

The  distance  between  the  side-walls  plays  a  significant  effect 
with  aspect  ratio  of  the  backward-facing  step  ducts.  De  Brederode 
and  Bradshaw  [12]  systematically  studied  the  side  wall  effects 
with  top  hat  free  stream  velocity  profile.  It  was  concluded  that  the 
three-dimensional  features  could  be  neglected  if  the  aspect  ratio  of 
the  step  was  larger  than  10.  Furthermore,  Hertzberg  and  Ho  [13] 
observed  that  the  presence  of  the  strong  spanwise  velocity  in  a 
small  aspect  ratio  rectangular  sudden  expansion  will  obviate  the 
formation  of  the  recirculation  zone  after  the  step  and  therefore  has 
a  zero  reattachment  length. 

Tylli  et  al.  [14]  studied  the  sidewall  effects  in  three-dimensional 
laminar  water  flow  over  a  backward-facing  step  experimen¬ 
tally  and  numerically.  Experiments  were  based  on  both  digital 
and  stereoscopic  particle  image  velocimetry;  a  spectral  element 
method  was  used  for  the  laminar,  transition,  and  turbulent  simula¬ 
tions.  The  geometry  has  an  expansion  ratio,  ER,  and  aspect  ratio,  AR, 
of  2  and  20,  respectively.  It  was  noticed  that  excellent  agreement 
between  the  3D  simulation  and  the  experiments  is  obtained  when 
Reynolds  number  reaches  Re ss  650.  Moreover,  it  was  found  that 
at  low  Reynolds  numbers  (Re  <400),  side  wall  effects  do  not  affect 
the  structure  of  laminar  flow  in  the  channel  mid  plane;  thus,  the 
mid  plane  field  is  accurately  predicted  by  two-dimensional  simu¬ 
lations.  In  addition,  it  was  found  that  at  higher  Reynolds  numbers, 
laminar  flow  is  characterized  by  sidewall  separation  and  the  for¬ 
mation  of  a  recirculation  zone,  which,  however,  does  not  penetrate 
up  to  the  channel  mid  plane.  The  wall-jet  intensity  and  flow  three- 
dimensionality  increase  with  Reynolds  number  for  laminar  flow, 
and  decrease  for  transitional  flow;  three-dimensional  variation  is 
hardly  noticed  in  the  time  average  field  of  turbulent  flow. 
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Armaly  et  al.  [15]  measured  the  velocity  for  three-dimensional 
laminar  separated  airflow  adjacent  to  a  backward-facing  step  using 
two-component  laser  Doppler  velocimeter.  The  backward-facing 
step,  with  a  height  of  s=  1.0  cm  has  an  aspect  ratio  of  AR  =  8,  and 
an  expansion  ratio  of  ER  =  2.02.  The  flow  measurements  were  cov¬ 
ered  Reynolds  number  range  between  98.5  <  Re  <525  to  develop 
a  clear  view  of  both  sidewall  and  the  step  effects,  and  compared 
with  the  previous  2D  studies.  It  was  determined  that  for  Re  <  98.5 
there  was  no  recirculation  flow  region  adjacent  to  the  sidewall  and 
at  Re<  190,  a  small  recirculation  flow  region  was  detected  in  the 
upper  corner  of  the  sidewall.  It  was  concluded  that  the  spanwise 
and  the  transverse  peak  velocities  move  to  the  center  of  the  duct 
as  the  distance  from  the  step  increases. 

Nie  and  Armaly  [16]  measured  the  velocity  adjacent  to  the 
bounding  walls  of  three-dimensional  backward-facing  step  flow 
using  Laser-Doppler  velocimeter.  The  study  was  performed  for  the 
purpose  of  mapping  the  boundaries  of  the  reverse  flow  regions  that 
develop  adjacent  to  the  walls  of  a  geometry  with  a  step  height 
s  =  1  cm,  aspect  ratio  AR  =  8,  and  an  expansion  ratio  ER  =  2.02  as  a 
function  of  Reynolds  number.  The  Reynolds  numbers  were  pre¬ 
sented  in  the  range  of  100  <  Re  <  8000  to  cover  the  three  types  of 
flow  regimes.  The  boundaries  of  the  reverse  flow  regions  were  iden¬ 
tified  by  locating  the  streamwise  coordinates  on  a  plane  adjacent 
to  the  bounding  walls  where  the  mean  streamwise  velocity  com¬ 
ponent  is  zero.  The  three-dimensional  reattachment  point  at  the 
stepped  wall  was  compared  with  the  two-dimensional  at  the  cen¬ 
ter  of  the  test  section.  It  was  noticed  that  the  three-dimension  flow 
results  are  slightly  higher  in  the  laminar  flow  regime  (Re  <400); 
significantly  lower  in  the  transition  flow  regime  (400  <  Re  <3400); 
and  slightly  lower  in  the  fully  turbulent  flow  regime  as  (Re  >  3400). 
It  was  noticed  that  the  size  of  the  reverse  flow  regions  increases  and 
moves  further  downstream  in  the  laminar  flow  regime;  decreases 
and  moves  upstream  in  the  transitional  flow  regime;  and  remains 
almost  constant  or  diminishes  in  the  turbulent  flow  regime;  as  the 
Reynolds  number  increases. 

The  literature  review  reveals  that  there  were  several  studies 
focused  on  the  heat  transfer  phenomena  over  the  backward-facing 
step  which  exhibit  an  influence  in  the  fluids  flow  behavior.  Aung 
[17],  and  Sparrow  and  Chuck  [18]  are  the  first  investigators  who 
established  and  reported  an  experimental  and  numerical  results 
on  heat  transfer  for  2D  laminar  air  flow  passing  a  backward-facing 
step  channel  heating  at  uniform  temperature  from  below.  In  their 
studies,  they  outlined  that  for  laminar  airflow  the  heat  transfer 
increases  monotonically  in  the  stream-wise  direction  and  quan¬ 
titatively  it  is  less  than  that  of  the  flat  plate  value.  It  was  concluded 
that  the  Nusselt  number  was  independent  of  the  Reynolds  number; 
the  local  Nusselt  number  distribution  begins  with  a  low  value  at 
the  step,  and  then  increases  monotonically  and  attains  a  maximum 
value  at  a  position  near  the  reattachment  point  then  it  decreases 
monotonically  towards  the  fully  developed  value. 

Khanafer  et  al.  [19]  analyzed  mixed  convection  of  laminar  pul¬ 
satile  flow  and  heat  transfer  past  a  backward-facing  step  in  a 
channel.  Fluid  flow  and  heat  transfer  characteristics  were  exam¬ 
ined  in  the  domain  of  the  Reynolds  number,  Richardson  number 
and  the  dimensionless  oscillation  frequency  as:  1 00  <  Re  <1000, 
1.78  x  10~3  <Ri<  10,  and  0.1<nr<5.  It  was  found  that  the 
elimination  of  the  separation  region  occurs  for  Ri>0.1,  and  as 
Reynolds  number  increases,  the  recirculation  zone  along  the 
heated  downstream  surface  becomes  larger.  The  thermal  bound¬ 
ary  layer  decreases  in  thickness  as  Reynolds  number  increases  and 
consequently  increases  the  heat  transfer  rate.  The  average  Nus¬ 
selt  number  and  the  length  of  the  recirculation  zone  decrease 
with  an  increase  in  Richardson  number.  It  was  also  noticed 
that  the  dimensionless  local  variation  of  skin  friction  coefficient 
along  the  heated  wall  increases  with  increasing  the  Richardson 
number. 


Chen  et  al.  [20]  simulated  numerically  the  turbulent  forced 
convective  flow  adjacent  to  a  two-dimensional  backward¬ 
facing  step.  The  objective  of  this  study  was  to  explore 
the  effects  of  step  height  on  turbulent  separated  flow  and 
heat  transfer.  The  stepped  wall  downstream  from  the  step 
was  set  to  be  uniform  and  constant  heat  flux  qw  =  270W/m2, 
while  other  walls  were  treated  as  adiabatic.  Reynolds  number 
and  duct’s  height  downstream  from  the  step  were  kept  constant 
at  Re0  =  28,000  and  H  =  0.19  m,  respectively.  On  the  other  hand, 
the  expansion  ratio  was  maintained  to  ER=1.11,  1.25,  and  1.67, 
respectively.  It  was  noticed  that  as  the  step  height  increases;  the 
primary  and  secondary  recirculation  regions,  magnitude  of  the 
maximum  turbulent  kinetic  energy  increase,  and  the  bulk  temper¬ 
ature  increases  rapidly  as  well.  Near  the  step  and  below  the  step 
height,  it  was  found  that  the  kinetic  energy  becomes  smaller  as  the 
step  height  increases.  Moreover,  it  was  found  that  the  maximum 
temperature  becomes  greater  as  the  step  height  increases.  The  peak 
values  of  the  transverse  velocity  component  become  smaller  as  the 
step  height  increases,  but  inside  the  recirculation  region;  the  skin 
friction  coefficient  becomes  less  significant  and  smaller  in  magni¬ 
tude  with  the  increase  of  the  step  height. 

The  effects  of  the  flow  and  heat  transfer  in  a  3D  backward-facing 
step  ducts  have  been  extensively  studied  using  different  parame¬ 
ters  by  Iwai  et  al.  [21],  Nie  and  Armaly  [22],  Armaly  et  al.  [23], 
Barbosa  Saldana  and  Anand  [24],  Barbosa  Saldana  et  al.  [5],  and 
Lan  et  al.  [25],  Iwai  et  al.  [21]  performed  a  numerical  simulation 
for  three-dimensional  mixed  convection  flows  over  a  backward¬ 
facing  step  at  low  Reynolds  number  in  order  to  investigate  the 
effects  of  the  duct  aspect  ratio.  Numerical  results  were  conducted 
for  an  aspect  ratio  range  between  4  <  AR  <  24,  and  Reynolds  num¬ 
ber  range  between  125<Re<375  to  study  its  effects  on  Nusselt 
number  and  the  skin  friction  coefficient  on  the  bottom  wall.  It  was 
found  that  an  aspect  ratio  of  16  at  least  is  needed  to  secure  2D 
region  in  the  central  part  of  the  duct  at  a  Reynolds  number  of  250. 
It  was  also  found  that  the  maximum  Nusselt  number  increases  with 
an  increase  of  the  aspect  ratio  and  Reynolds  number.  The  skin  fric¬ 
tion  coefficient  is  found  to  be  an  analogous  to  the  case  of  Nusselt 
number  distribution. 

Nie  and  Armaly  [22]  presented  a  numerical  analysis  for  a 
three-dimensional  incompressible  laminar  forced  convection  flow 
adjacent  to  backward-facing  step  in  rectangular  duct.  The  effects 
of  step  height  on  the  flow  and  heat  transfer  characteristics  were 
the  main  objectives  of  this  study.  The  geometry,  thermophysical 
properties,  and  flow  condition  are  based  on  the  available  measure¬ 
ments  of  Armaly  et  al.  [23],  The  simulations  were  performed  at 
Re  =  343  for  different  step  heights  (s  =  0.008,  0.010,  and  0.012  m), 
while  the  other  geometries  were  kept  constant.  It  was  found  that 
increasing  the  step  height  increases  the  reattachment  length,  the 
Nusselt  number,  the  size  of  the  sidewall  reverse  flow  region,  and 
the  general  three-dimensional  features  of  the  flow.  In  addition, 
as  the  step  height  increases,  the  distance  between  the  step  and 
the  point  where  the  impingement  flow  reaches  the  stepped  wall 
downstream  increases,  and  a  secondary  recirculation  flow  region 
develops  adjacent  to  the  bottom  corner  of  the  step.  For  the  case 
of  small  step  height  (s  =  0.008  m),  the  minimum  Nusselt  number  is 
located  near  the  bottom  corner  between  the  step  and  the  center- 
line  of  the  duct.  As  the  step  height  increases,  the  minimum  Nusselt 
number  is  appeared  to  move  towards  near  the  bottom  corner  of 
the  step  and  the  sidewall.  It  was  found  that  the  friction  coefficient 
increases  in  the  streamwise  direction  along  the  centerline  of  the 
duct  inside  the  primary  recirculation  flow  region,  but  it  decreases 
outside  the  recirculation  region  with  increasing  step  height. 

Barbosa  Saldana  and  Anand  [24]  studied  numerically  forced  con¬ 
vective  airflow  over  a  three-dimensional  backward-facing  step.  The 
horizontal  backward-facing  step  channel  was  set  to  have  a  step 
height  s=  1  cm,  an  aspect  ratio  of  8  and  an  expansion  ratio  of  2.  The 
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stepped  wall  was  fixed  to  a  constant  heat  flux  gw  =  50  W/m2  while 
the  other  walls  were  considered  as  insulated.  The  Reynolds  num¬ 
bers  were  set  to  be  in  the  range  of  98.5  and  512.  The  local  Nusselt 
number  along  the  bottom  wall  was  found  to  lie  in  the  vicinity  of  the 
xu-line  and  the  xw-line  points  of  intersection  where  the  shear  stress 
along  the  bottom  wall  is  equal  to  zero.  Furthermore,  the  velocity 
profiles  revealed  that  for  Reynolds  number  greater  than  343  the 
flow  does  not  reach  fully  developed  conditions  at  the  channel  exit. 

Barbosa  Saldana  et  al.  [5]  simulated  numerically  laminar  mixed 
convective  flow  over  a  three-dimensional  horizontal  backward¬ 
facing  step.  All  the  walls  are  treated  adiabatically  except  the  bottom 
wall  of  the  channel  is  subjected  to  a  constant  high  temperature. 
An  airflow  of  fixed  Reynolds  number  Re  =  200  enters  hydrodynam- 
ically  fully  developed  and  isothermally  in  a  duct  with  an  aspect 
ratio  of  4,  expansion  ratio  of  2,  and  total  length  of  52  times  the 
step  height  L  =  52S.  It  was  reported  that  as  Richardson  number  (Ri) 
increases;  the  location  of  xu-line  displaces  to  the  upstream  (higher 
velocity  components),  the  recirculation  zone  is  shortened,  verti¬ 
cal  size  of  the  recirculation  zone  reduced,  and  negative  values  of  u 
velocity  ( streamwise  direction)  next  to  the  bottom  wall  of  the  chan¬ 
nel  increases.  In  addition,  as  the  Richardson  number  increased;  the 
location  of  the  maximum  averaged  Nusselt  number  distribution 
was  moved  and  the  maximum  averaged  streamwise  shear  stress 
distribution  moved  its  location  further  upstream. 

Lan  et  al.  [25]  simulated  the  turbulent  forced  convection  adja¬ 
cent  to  three-dimensional  backward-facing  step  in  a  rectangular 
duct  using  a  ( k-e  -  £-f)  turbulent  model.  The  aim  of  the  study 
was  to  present  a  clear,  realistic,  and  more  accurate  investigation 
near  the  walls  to  predict  the  turbulent  heat  transfer  separated 
and  wall-bounded  flows  better  than  the  commonly  used  two- 
equation  turbulence  model.  A  User  Defined  Function  (UDF)  was 
implemented  using  FLUENT-CFD,  and  two-dimensional  benchmark 
problems  were  simulated  to  validate  the  three-dimensional  simu¬ 
lation  code.  The  geometry  has  an  expansion  ratio  and  a  step  height 
of  1 .48  and  4.8  mm,  respectively.  The  three  aspect  ratios  of  3,  8  and 
infinity  (2-D  simulation)  were  considered  for  studying  its  effect  on 
the  flow  and  heat  transfer.  The  effect  of  Reynolds  number  in  the 
range  of  20,000  <  Re  <50,000  was  also  considered.  The  effects  of 
the  Reynolds  number  and  aspect  ratio  on  the  flow  reattachment 
was  minimal  in  the  range  of  parameters  that  were  examined,  but 
they  were  shown  to  have  significant  influence  on  the  heat  trans¬ 
fer.  The  reattachment  length  was  found  to  be  constant  in  central 
portion  of  duct  width,  but  starts  to  increase  as  it  approaches  the 
side  walls  and  reaches  its  maximum  value  on  the  side  walls.  It  was 
noticed  in  the  primary  recirculation  flow  region  and  in  the  vicin¬ 
ity  of  the  recirculation  line  a  significant  spanwise  flow  developed. 
The  maximum  local  Nusselt  number  was  found  to  be  developed  in 
the  region  where  the  spanwise  velocity  is  maximum  adjacent  to 
the  heated  wall  close  to  the  sidewall.  The  bulk  air  temperature  was 
found  to  decrease  as  the  Reynolds  number  increases,  but  it  was  not 
influenced  significantly  by  the  changes  in  the  aspect  ratio. 

3.2.  Inclined  flow 

The  backward-facing  step  channel  orientation  effects  attracted 
some  researchers  such  as  Lin  et  al.  [26],  Lin  et  al.  [27],  Flong  et  al. 
[28],  Abu-Mulaweh  et  al.  [29],  and  Iwai  et  al.  [30]  to  investigate 
the  separation  region  and  the  reattachment  point  with  respect 
to  other  parameters.  The  literature  declared  that  the  first  effort 
dealt  with  backward-facing  step  in  inclined  duct  was  introduced 
by  Lin  et  al.  [26],  They  examined  numerically  the  effects  of  inclina¬ 
tion  angle  (0°  <  cp<  90°  from  the  vertical)  on  the  two-dimensional 
mixed  convective  laminar  flow  and  heat  transfer  in  a  duct  with 
a  backward-facing  step  under  buoyancy-assisting  flow  conditions. 
The  case  study  was  taken  from  a  previous  study  done  by  the  same 


author  for  a  vertical  duct  with  a  backward-facing  step  Lin  et  al. 
[27],  The  geometry  has  an  expansion  ratio  and  step  height,  of  2  and 
0.48  cm,  respectively.  The  straight  wall  at  the  top  of  the  duct  was 
maintained  at  the  inlet  air  temperature  of  20  °C,  the  downstream 
stepped  wall  was  maintained  to  a  uniform  temperature  of  35  °C, 
while  the  upstream  wall  and  the  step  were  maintained  adiabatic. 
The  Reynolds  number  was  maintained  at  50.  It  was  found  that  the 
recirculation  size  increases  as  the  orientation  transfer  from  the  ver¬ 
tical  axis.  The  streamwise  component  of  the  Grashof  number  has 
a  strong  influence  on  the  flow,  but  the  transverse  component  has 
a  small  effect  on  the  velocity  and  temperature  profiles  in  the  fully 
developed  regime.  The  substantial  fact  elaborates  that  the  reattach¬ 
ment  length,  fluid  temperature,  and  the  shear  stress  at  the  unheated 
wall  increase  as  the  inclination  angle  increases  from  the  vertical.  In 
contrary,  for  Nusselt  number  for  both  heated  and  unheated  walls, 
the  streamwise  buoyancy  force,  and  the  shear  stress  at  the  heated 
wall  are  decreased. 

Hong  et  al.  [28]  were  the  first  who  investigated  the  buoyancy 
opposing  effects  in  an  inclined  channel  with  a  backward-facing 
step.  They  studied  numerically  the  influence  of  the  relevant  inclina¬ 
tion  angle  and  Prandtl  number  for  two-dimensional  laminar  mixed 
convective  airflow  for  buoyancy  assisting  and  opposing  flow  condi¬ 
tions.  The  geometry  has  an  expansion  ratio  and  step  height,  of  2  and 
0.0048  m,  respectively.  The  straight  wall  at  the  top  of  the  duct  was 
maintained  at  the  inlet  air  temperature  of  293  K,  the  downstream 
stepped  wall  was  maintained  to  a  uniform  heat  flux  of  200  W/m2, 
while  the  upstream  wall  and  the  step  were  maintained  adiabatic. 
The  Reynolds  number  was  fixed  at  100,  and  the  Grashof  number 
was  fixed  at  609.  The  inclination  angle  was  varied  from  0°  to  360°, 
while  other  parameters  were  kept  constant.  It  was  found  that  as  the 
inclination  angle  increases  from  0°  to  1 80°,  the  reattachment  length 
increases,  but  the  wall  friction  coefficient  and  the  Nusselt  number 
at  the  heated  wall  decrease.  Moreover,  the  transverse  buoyancy 
force  was  noticed  to  have  a  less  significant  effect  on  the  flow  and 
heat  transfer  downstream  the  reattachment  point,  while  it  has  a 
strong  effect  in  the  recirculation  region.  The  Prandtl  number  was 
varied  in  the  range  of  0.07  <Pr<  100,  while  the  inclination  angle 
was  fixed  at  0°.  It  was  found  that  as  the  Prandtl  number  increases, 
the  Nusselt  number  and  the  reattachment  length  increase,  but  the 
wall  friction  coefficient  decreases. 

The  understanding  of  the  separated  zone  and  reattachment  pro¬ 
cess  was  poor  for  the  inclined  channel  with  a  backward-facing 
step  due  to  the  inability  to  predict  simple  reattaching  flows  over 
a  wide  range  of  parameters.  The  first  experimental  effort  dealt 
with  backward-facing  step  mounted  in  inclined  channels  is  pre¬ 
sented  by  Abu-Mulaweh  et  al.  [29].  They  measured  experimentally 
and  simulated  numerically  the  influence  of  buoyancy-assisting 
laminar  mixed  convection  in  boundary-layer  flow  over  horizontal 
and  inclined  two-dimensional  backward-facing  step.  The  objec¬ 
tives  of  the  study  were  to  determine  the  reattachment  length 
and  the  onset  of  vortex  instability  for  different  wall  temperatures 
(0  <  AT <  30  °C),  free  stream  velocities  (0.285  <  Uoo  <  0.7  m/s),  step 
heights  (0.35  < s  <  0.8  cm),  and  inclination  angles  (30°  <cp<  90°)  as 
measured  from  the  vertical  direction.  The  geometry  of  the  test 
section  was  consisted  of  a  heated  constant  temperature  down¬ 
stream  section,  adiabatic  upstream  section  and  step,  and  the  top 
section  was  a  free  stream.  It  was  noticed  in  horizontal  case  that 
the  buoyancy  force  has  a  negligible  effect  on  the  velocity  and  tem¬ 
perature  distributions,  but  it  significantly  influences  the  onset  of 
instability.  On  the  other  hand,  it  was  found  that  as  the  inclina¬ 
tion  angle  increases  from  the  vertical,  the  streamwise  component 
of  the  buoyancy  force  decreases,  which  causes  an  increase  in 
the  velocity  gradient  and  a  decrease  in  the  temperature  gradient 
at  the  wall  in  the  recirculation  region.  It  was  concluded  that  as 
the  inclination  angle  increases  from  the  vertical,  the  local  Nus¬ 
selt  number  decreases,  while  the  reattachment  length  and  the 
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location  of  the  maximum  Nusselt  number  behind  the  backward¬ 
facing  step  increase.  It  was  noticed  that  the  maximum  Nusselt 
number  occurs  downstream  of  the  reattachment  point.  Moreover, 
the  space  between  the  maximum  Nusselt  number  location  and  the 
reattachment  decreases  slowly  as  the  inclination  angle  increases. 

As  it  has  been  introduced  in  the  previous  sections,  all  of  the  exist¬ 
ing  researches  for  the  inclined  channel  with  a  backward-facing  step 
are  focused  on  two-dimensional  situations,  though  most  of  the  real¬ 
istic  applications  are  practically  dominated  by  three  dimensional 
geometries.  It  is  declared  from  the  literature  that  the  first  three- 
dimensional  study  is  presented  by  Iwai  et  al.  [30].  They  simulated 
numerically  the  effect  of  inclination  angle  in  three-dimensional 
mixed  convective  flow  over  a  backward-facing  step  in  a  rectan¬ 
gular  duct.  The  downstream  wall  was  kept  with  a  constant  heat 
flux,  while  other  walls  were  kept  adiabatic.  Parameters  such  as 
Reynolds  number,  expansion  ratio,  and  aspect  ratio  were  kept  con¬ 
stant,  at  125,  2,  and  16,  respectively.  The  effects  of  two  inclination 
angles;  <f>\  pitch  angle  changing  in  the  range  of  0°  <  0i  <  360°,  and 
<p2  rolling  angle  changing  in  the  range  of  0°<(/>2<180°  are  the 
main  interest  of  this  study.  The  three  positions  of  the  reattachment 
point  X, ,  the  peak  Nusselt  number  point  Xp  and  the  downstream 
end  of  the  secondary  recirculation  region  Xs  on  the  centerline  of 
the  heated  wall  are  shifted  with  tf>\.  It  was  found  that  when  (/q 
is  changed  while  keeping  <f>2  to  be  zero  degree,  positions  of  the 
maximum  Nusselt  numbers  are  symmetrically  obtained  at  the  posi¬ 
tions  on  the  heated  wall  near  the  two  side  walls,  similar  to  the 
cases  of  pure  forced  convection.  The  maximum  Nusselt  number 
was  appeared  at  the  most  upstream  position  in  the  case  of  </q  =  0°, 
and  the  minimum  at  the  most  downstream  position  in  the  case  of 
<pl  =  180°.  On theother  hand,  itwas  noticed  thatwhen(/>2  is  changed 
while  keeping  4> i  to  be  90°,  the  flow  and  thermal  fields  become 
asymmetric  about  the  duct  centerline.  It  was  observed  that  4>2 
has  little  effect  on  the  spanwise  position  of  the  maximum  Nusselt 
number. 

3.3.  Vertical  flow 

It  is  indicated  from  the  literature  that  the  first  effort  of  study¬ 
ing  the  effects  of  vertical  backward-facing  step  was  introduced 
by  Lin  et  al.  [31].  They  examined  and  reported  the  numerical 
results  of  laminar  mixed  convection  airflow  and  heat  transfer, 
buoyancy-assisting  through  a  two-dimensional  vertical  duct  with 
a  backward-facing  step.  The  study  was  conducted  to  cover  the 
domain  from  pure  forced  convective  flow  to  the  inlet  starved 
convective  flow  where  the  average  inlet  velocity  is  smaller  than 
the  corresponding  natural  convective  value.  The  geometry  dimen¬ 
sions  and  boundary  conditions  were  based  on  previous  work  done 
by  Armaly  et  al.  [7[.  The  influence  of  the  temperature  difference 
between  the  heated  wall  downstream  of  the  step  and  free  stream 
in  the  range  of  1  <  AT <  75  °C  was  detected.  It  was  found  that  as  the 
wall  temperature  increases,  the  reattachment  length  decreases  and 
the  size  of  the  secondary,  inner  recirculation  region  that  develops  at 
the  corner  between  the  heated  wall  and  the  step  increases.  On  the 
other  hand,  the  buoyancy  level  decreases  for  a  fixed  wall  tempera¬ 
ture  as  the  Reynolds  number  increases.  It  was  noticed  that  the  force 
of  the  pressure  drop  Fp  is  lower  than  the  force  of  the  kinetic  energy 
Fk,  where  a  reverse  flow  could  be  developed  from  the  downstream 
section  of  the  duct  along  the  cooler  wall  as  outlined  by  Aung  and 
Worku  [32].  However,  it  was  found  that  the  disappearance  of  the 
reattachment  from  the  heated  wall  will  occur  only  in  the  starved 
flow  regime.  The  Nusselt  number  was  found  to  be  increased  as  the 
wall  temperature  increases  and  the  Reynolds  number  decreases. 
Both  the  peak  Nusselt  number  length  X„,  and  the  reattachment 
length  of  the  main  recirculation  region,  Xr,  decrease,  while  the 
length  of  the  recirculation  region,  X0,  increases  as  the  value  of 
the  buoyancy  parameter  increases.  Moreover,  the  peak  Nusselt 


number  occurs  downstream  of  the  reattachment  point  and  the 
spacing  increases  as  the  buoyancy  parameter  increases. 

The  effects  of  step  height  in  inclined  ducts  that  presented  before 
in  this  literature  did  not  include  the  vertical  orientation  (<p  =  0°). 
This  gap  of  research  motivated  Hong  et  al.  [4,33]  to  investigate 
the  effects  of  other  parameters  due  to  different  step  height.  They 
numerically  studied  the  effects  of  the  mixed  convection  flow  and 
heat  transfer  through  a  two-dimensional  vertical  duct  having  a 
backward-facing  step.  The  downstream  wall  from  the  step  was 
subjected  to  a  uniform  heat  flux  (UHF),  while  the  upstream  wall 
and  the  step  by  itself  were  adiabatic,  and  the  straight  top  wall  was 
maintained  at  a  uniform  temperature  that  is  equal  to  the  inlet  fluid 
temperature.  Furthermore,  they  examined  the  effects  of  different 
ranges  of  parameters  such  as  Grashof  number,  0  <  Crs*  <  24,000, 
expansion  ratio,  1.25<£R<4,  and  Reynolds  number,  0<Re<150. 
It  was  found  that  Grashof  number  has  a  significant  influence  on 
the  flow  and  heat  transfer;  when  Grashof  number  increases,  the 
wall  friction  coefficient  and  the  local  Nusselt  number  on  the  down¬ 
stream  heated  wall  increase,  and  verse  visa  for  the  reattachment 
length.  It  was  noticed  that  the  reattachment  length  increases  as 
the  expansion  ratio  increases  for  less  than  2.25,  but  for  values 
greater  than  2.25  the  reattachment  length  decreases.  Moreover,  it 
was  found  that  increasing  the  expansion  ratio  increases  the  fully 
developed  value  of  the  wall  friction  coefficient  and  the  distance  for 
the  flow  to  reach  the  fully  developed  regime,  but  it  decreases  the 
peak  and  the  fully  developed  Nusselt  numbers.  They  have  devel¬ 
oped  correlation  equations  for  reattachment  length  and  the  peak 
Nusselt  number  as  shown  in  Table  1. 

Baek  et  al.  [34]  studied  experimentally  and  numerically  the 
buoyancy-assisting  laminar  mixed  convection  airflow  and  heat 
transfer  over  a  vertical,  two-dimensional  backward-facing  step. 
Different  free-stream  velocities,  wall  temperature  differences,  and 
step  heights  in  the  range  of  (0.37  <  uo<0.72  m/s),  (10  <  AT<30°C), 
and  (0.38  <s<  1  cm),  respectively,  were  considered  to  visualize 
the  flow  and  heat  transfer  behavior.  It  was  found  that  for  a  fixed 
inlet  velocity,  the  reattachment  length,  Xr,  decreases  as  the  wall 
temperature  increases.  Similarly,  for  a  fixed  wall  temperature,  the 
reattachment  length,  Xr,  and  the  volume  of  the  recirculation  region 
decrease  as  the  inlet  velocity  decreases.  It  was  noticed  that  the 
recirculation  region  disappears  and  detaches  from  the  heated  wall 
and  attaches  only  to  the  adiabatic  step.  Moreover,  for  the  heat 
transfer  visualization,  it  was  found  that  the  heat  transfer  rate  is 
seen  to  increase  as  the  temperature  difference  increases.  The  fluid 
temperature  in  the  recirculation  region  behind  the  step  (X<Xr) 
was  considerably  higher  than  the  temperature  downstream  of  the 
reattachment  point  (X>Xr).  As  the  buoyancy  force  increases,  the 
magnitude  of  the  Nusselt  number  increases  and  the  location  of  its 
peak  moves  closer  to  the  step.  An  effort  was  made  by  the  authors 
to  modify  the  existing  empirical  correlation  for  the  reattachment 
length  in  the  laminar  forced  convection  flow,  as  given  by  Goldstein 
et  al.  [6],  to  reflect  also  the  influence  of  step  height  and  buoyancy 
parameter.  It  was  found  that  the  location  of  the  maximum  Nus¬ 
selt  number,  Xn,  and  the  reattachment  length,  Xr,  decreases  as  the 
buoyancy  force  increases. 

Abu-Mulaweh  et  al.  [4,35]  examined  the  influence  of  buoyancy- 
assisting  laminar  flow  and  heat  transfer  over  a  2D  vertical 
backward-facing  step.  The  boundary  conditions  and  general  step 
geometry  were  taken  from  Baek  et  al.  [34]  except  the  downstream 
heated  wall  was  maintained  to  a  uniform  heat  flux  in  the  range 
between  102<qw<290W/m2.  As  the  buoyancy  force  increases, 
the  reattachment  length  Xr  decreases  but  the  local  Nusselt  num¬ 
ber  increases  and  the  location  of  its  maximum  values  moves  closer 
to  the  step.  Moreover,  it  was  found  that  the  local  Nusselt  number  to 
be  lower  in  the  region  where  X<Xr  and  higher  in  the  region  where 
X>Xr  for  higher  free  stream  velocity.  Correlation  equations  for  the 
reattachment  length  were  reported  in  Table  1. 


Table  1 

Summary  of  the  research  done  on  backward  facing  step. 


Orientation  Heating  Fluid  Dim  Study  Ranges  Xr  correlation  equations  Xo  correlation  equations  Xn  correlationequations  NuSiTnax  correlation 

type  study  equations 


Denham  and 
Patrick  [7] 

Hor 

Water 

2D 

E 

50  <  Re  <250 

ER  =  3 

AR^20 

Armaly  et  al.  [8] 

Hor 

Air 

2D 

E&N 

70  <  Re  <8000 
ER=  1.9423 

AR  =  36 

Shih  and  Ho  [10] 

Hor 

- 

Water 

2D&3D 

E 

AR  =  3 

Laminar 

- 

- 

- 

- 

Chiang  et  al.  [11] 

Hor 

Oil 

3D 

E 

Re  <  800 

ER=  1.9432 

AR  =  10.6022 

Tylli  et  al.  [14] 

Hor 

Water 

3D 

E&N 

Lam<Re<Tur 

ER  =  2 

AR  =  20 

Armaly  et  al.  [15] 

Hor 

Air 

3D 

E&N 

98.5  <  Re  <525 

ER  =  2.02 

AR  =  8 

Nie  and  Armaly 
[16] 

Hor 

Air 

3D 

E 

1 00  <  Re  <8000 
ER=2.02 

AR  =  8 

Khanafer  et  al. 

[19] 

Hor 

UWTa 

Air 

2D 

Na 

100  <Re<  1000 
1.78E-3  <Rz  <10 
0.1  <UT  <  0.5 

AR  =  8 

Chen  et  al.  [20] 

Hor 

UHFa 

Air 

2D 

N 

Re  =  28000 

1.11  <ER<1.67 

- 

- 

- 

- 

Abu-Nada [169] 

Hor 

UWT 

NFa 

2D 

N 

200  <  Re  <600 

0<  UT<  0.2 

ER  =  2 

Iwai  et  al.  [21] 

Hor 

UWT 

Air 

3D 

N 

125  <  Re  <375 

4 

<  AR  <  24 

Nie  and  Armaly 
[22] 

Hor 

UHF 

Air 

3D 

N 

Re =343 

8  <  s  <  1 2  mm 

- 

- 

- 

- 

Saldana  and 

Anand  [24] 

Hor 

UHF 

Air 

3D 

N 

98.5  <Re<  512 

ER  =  2 

AR  =  8 

Saldana  et  al.  [5] 

Hor 

UWT 

Air 

3D 

N 

Re  =  200 

ER  =  2 

AR  =  4 

Lan  et  al.  [25] 

Hor 

UHF 

Air 

3D 

N 

2E4  <Re<  5E4 

ER=  1.48 

3  <  AR  <  oo 

Lin  et  al.  [26] 

Inca 

UWT 

Air 

2D 

N 

0 °  <  0  <  90° 

Re  =  50 

ER  =  2 

Hong  et  al.  [28] 

Inc 

UHF 

Air 

2D 

N 

0°  <0<36O° 

0.07  <Pr<  100 

Re =100 

Crs*  =  609 

ER  =  2 

Xr  =  [5.462 -1.62  cos(0)] 
x  [0.95  + 0.05  cos(20)][l 
+  O.O2sin(0)] 

Xr  =  4.9 -1.7 
exp(-Pr/0.4) 

X0  =  MAX{0.4462  cos(0)[l 
-O.5sin(0)],O} 

X0  =  0.336- 0.001 65Pr 

Xn  =  [5.472  -  0.6535  cos(0)][O98 
+  0.02  cos(20)][l  +O.O2sin(0)] 

Xn  =  3.53  +  1. 7  IPr, 

0.07  <Pr<  0.712 

Xn  =  4.753-  0.009 IPr2, 

0.712  <Pr<  100 

Nus>max=[l.Ol-O.Olcos(20)] 
[1.752  +  0.151  cos(0)] 

NuSlmax  =  [0.52  +  1.29VPr 
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Abu-Mulaweh 
et  al.  [29] 

Inc 

UWT 

Air 

2D 

E&N 

30°  <  0  <  90° 

50<Res  <300 

3.5  <  s  <  8  mm 

0<  AT<30°C 

Xr  =  0.023(8s/s) 

Res  +f[s/Xi)  =  0.023 
[5  (xi/s)Rexr1/2]Res+./[s/x,) 
j[s/Xi)  =  8995.8  (s/Xj)2  -  114.2 
(s/x  0  — 3.4 

Iwai  et  al.  [30] 

Inc 

UHF 

Air 

3D 

N 

0°  <  0i  <  360° 

0°  <  02  <  1 80° 

Re= 125 

ER  =  2,  AR  =  16  Ri  =  0.03 

Lin  et  al.  [31] 

Vera 

UWT 

Air 

2D 

N 

0  =  0° 

Res  =  50, 100 

1  <  AT<75°C 

ER  =  2 

Hong  et  al.  [33] 

Ver 

UHF 

Air 

2D 

N 

0  =  0° 

0  <  Res  <150 

0  <  Grs*  <  2400 

1.25  <ER<4 

Xr  =xr , forced  {1  +(7.188  —4.0471?  - 

-1.758Res)j} 

.forced  =  3o(R)+  ai  (R)Res 
a0(R)  =  -  1 .252  + 1 .3451?  -  0.1 681R2 
a,  (R)  =  - 0.02751 

0.04572R-  0.00641 4i?2 

NuSiniax  =  1 .658  +  0.0001 483Grs 
Nus,max  =  4.862- 2.1 72R 
+  0.3292R2 

Baek  et  al.  [34] 

Ver 

UWT 

Air 

2D 

E&N 

0  =  0° 

0.37  <  Uoo  <  0.72 

10<  AT<30°C 

3.8  <s<  10mm 

Xr  =  (2.24  +  0.022Res) 
exp[-42  ^s/x,)-0186] 

Abu-Mulaweh 
et  al.  [35] 

Ver 

UHF 

Air 

2D 

E 

0  =  0° 

50<Res  <330 

0  <  £  <  1 .4E-2 

3.5  <  s  <  8  mm 

Xr=Xriforcedexp(-7.5^/2)  # 

Xr/orc£  =  0.023(8s/s)Res+Ks/Xi) 

Ks/x,)  =  8995.8/[s/x,)2  - 114.2 
(s/x,)- 3.4 

# 

# 

Abu-Mulaweh 
et  al.  [36] 

Ver 

UWT 

Air 

2D 

E&N 

0=180° 

0.41  <  Uoo  <  0.66 

0<  AT<30°C 
4<s<8mm 

Iwai  et  al.  [38] 

Ver 

UHF 

Air 

3D 

N 

0  =  0 

0  <  ^  <  0.12 

Re= 125 

ER  =  2,  AR  =  16 

Abu-Mulaweh 
et  al.  [39] 

Ver 

UWT 

Air 

2D 

E 

0  =  0 

0  <  Uoo  <  0.41 

5  =  22  mm 

Abu-Mulaweh 
et  al.  [40] 

Ver 

UWT 

Air 

3D 

E 

0  =  0 

Uoo  =  0.41  m/s 

0  <  s  <  22  mm 

a  Hor,  horizontal;  Inc,  Inclined;  Ver,  vertical;  UHF,  uniform  heat  flux;  UWT,  uniform  wall  temperature;  NF,  nanofluids;  E,  experimental;  N,  numerical;  #,  not  identified;  no  correlation  equation. 
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As  introduced  in  the  previous  section  in  this  literature,  the 
researchers  were  concentrated  on  buoyancy-assisting  vertical 
flow.  The  first  effort  to  study  the  influence  of  the  buoyancy- 
opposing  flow  was  conducted  by  Abu-Mulaweh  et  al.  [36].  They 
performed  experimental  and  numerical  studies  of  laminar  flow 
over  a  vertical  backward-facing  step  to  measure  the  effect  of  the 
buoyancy-opposing  on  the  airflow  and  heat  transfer  behavior.  The 
experimental  setup,  test  section  and  backward-facing  step  geom¬ 
etry,  and  the  boundary  conditions,  are  identical  to  that  used  by 
Abu-Mulaweh  et  al.  [35]  except  that  the  orientation  of  the  air  tunnel 
is  changed  by  180°  to  create  buoyancy-opposing  flow  conditions. 
It  was  noticed  for  low  buoyancy  levels  Ri  <  4.4  x  1 0-3 ;  the  length  of 
the  recirculation  region  downstream  of  the  backward-facing  step 
increases  rapidly  as  the  buoyancy  level  increases.  However,  for 
low  buoyancy  levels  above  4.4  x  10-3;  the  flow  inside  the  recir¬ 
culation  region  remains  laminar  but  the  length  of  the  recirculation 
region  decreases  rapidly  due  to  the  effect  of  the  transition  of  the 
flow  from  laminar  to  turbulent  and  the  flow  becomes  turbulent 
downstream  of  the  recirculation  region.  The  effect  of  temperature 
difference  between  the  heated  wall  downstream  the  step  and  the 
free  stream  in  the  range  of  0  <  AT <  30  °C  was  studied.  It  was  found 
that  as  the  temperature  difference  AT  increases,  the  reattachment 
length  increases  and  the  Nusselt  number  decreases. 

Abu-Mulaweh  et  al.  [4,37]  studied  the  effects  of  heating  the 
upstream  wall  and  the  step  of  two-dimensional,  vertical  backward¬ 
facing  step  geometry  on  the  flow  and  heat  transfer  characteristics 
in  the  region  downstream  of  the  step.  Four  different  upstream  wall 
and  step  heating  conditions  were  examined.  The  first  case  Tem¬ 
perature  Boundary  Condition  (TBC)  was  TBC-1,  which  corresponds 
to  /wi  —  Too  =  Tw  —  Too,  the  second  case  was  TBC-2  corresponds  to 
Twl  -Too  =  2(Tw-Too),  the  third  case  was  TBC-1/2  corresponds  to 
TW1  -  Too  =  (Tw  -  Too)/2,  and  the  fourth  case  was  TBC-A  corresponds 
to  adiabatic  upstream  wall  and  step  (where  Tw  is  the  downstream 
wall  temperature,  Twl  is  the  upstream  wall  and  step  temperature, 
and  Too  is  the  free  stream  temperature).  It  was  found  that  the  recir¬ 
culation  region  remains  attached  to  the  heated  wall  downstream 
of  the  step  when  the  upstream  wall  and  the  step  are  heated  to 
the  same  or  different  uniform  temperatures  as  the  downstream 
wall.  The  detach  of  the  recirculation  region  from  the  heated  wall 
downstream  of  the  step  may  occur  at  high  buoyancy  levels  when 
the  upstream  wall  and  step  are  maintained  adiabatic.  For  higher 
upstream  heating  levels,  it  was  found  that  the  reattachment  length 
decreases  as  the  buoyancy  level  increases.  Furthermore,  it  was 
noticed  that  when  the  upstream  wall  and  the  step  are  heated  to 
a  higher  temperature  than  the  downstream  wall,  the  velocity  dis¬ 
tribution  downstream  of  the  step  may  exhibit  an  overshoot  over  its 
free  stream  value.  It  was  also  noticed  that  the  heat  transfer  from 
the  downstream  wall  to  the  fluid  decreases  as  the  upstream  wall 
and  the  step  are  heated. 

Iwai  et  al.  [38]  presented  a  numerical  study  to  investigate  the 
effects  of  buoyancy-assisting  laminar  mixed  convective  flow  and 
thermal  fields  through  a  3D  vertical  duct  having  backward-facing 
step.  The  downstream  wall  was  maintained  to  a  uniform  heat  flux, 
while  the  rest  of  the  walls  were  assumed  to  be  adiabatic.  The  simu¬ 
lations  were  carried  out  for  buoyancy  parameter  range  between 
0<?<0  .12,  other  parameters  such  as  Reynolds  number,  aspect 
ratio,  expansion  ratio,  were  maintained  to  be  fixed  at  125,  16,  2, 
respectively.  Iwai  and  his  coworkers  have  noticed  that  2D  are  effec¬ 
tive  in  predicting  such  type  of  phenomena  for  buoyancy  parameter 
f  <  0.06,  for  values  greater  than  this,  a  3D  computational  study 
becomes  necessary.  The  reattachment  point  and  the  peak  Nusselt 
number  move  closer  to  the  step  as  the  buoyancy  force  increases. 
Moreover,  the  secondary  recirculation  region  which  developed  at 
the  corner  between  the  step  and  the  stepped  downstream  wall 
increases  in  size  as  the  buoyancy  force  increases.  It  was  found  that 
the  maximum  Nusselt  numbers  are  located  symmetrically  near  the 


side  walls  far  from  the  center  plane  of  the  downstream  wall.  In 
addition,  the  maximum  Nusselt  number  was  noticed  to  be  down¬ 
stream  of  the  reattachment  point  and  the  distance  between  them 
increases  with  increasing  buoyancy  force. 

The  research  to  study  effects  of  separation  and  reattachment 
phenomenon  does  not  concern  on  laminar  flow.  Abu-Mulaweh 
et  al.  [39]  examined  the  effect  of  turbulent  mixed  convection  over 
a  2D  vertical  backward-facing  a  step.  The  step  geometry  consists 
of  an  adiabatic  backward-facing  step,  an  upstream  wall  and  a 
downstream  wall  heated  at  constant  and  uniform  temperatures. 
A  free  stream  air  was  used  as  a  medium  in  the  velocities  range  of 
0  <  Uoo  <  0.41  m/s  for  a  step  height  of  22  mm,  and  a  temperature  dif¬ 
ference  (AT=30°C)  between  the  heated  walls  and  the  free  stream. 
The  measurements  for  different  free  stream  velocities  (Uoo  =  0, 0.23, 
0.41  m/s)  were  carried  out  at  one  streamwise  location  upstream  of 
the  step  (x  =  -5  cm)  and  four  different  streamwise  locations  down¬ 
stream  of  the  step  (x  =  3.5,  6,  10,  25  cm).  It  was  indicated  that  the 
streamwise  velocity  increases  as  the  free  stream  increases  and 
vice  versa  for  the  temperature  distribution  for  both  upstream  and 
downstream  locations.  The  reattachment  length  was  noticed  to  be 
increased  with  increasing  the  free  stream  velocity.  The  magnitudes 
for  the  intensities  of  streamwise  and  transverse  velocity  fluctua¬ 
tions  were  found  to  be  larger  for  natural  convection  than  for  mixed 
convection  on  the  upstream  location.  While,  the  intensities  tem¬ 
perature  fluctuations  decrease  rapidly  for  mixed  convective  flow. 
On  the  other  hand,  the  magnitudes  of  the  intensities  of  both  veloc¬ 
ity  and  temperature  fluctuations  are  smaller  for  mixed  convection 
than  natural  convection  at  the  four  downstream  locations.  It  was 
found  that  the  reattachment  length  increases  while  the  Nusselt 
number  from  the  downstream  heated  wall  decreases  as  the  free 
stream  increases,  and  as  a  result  the  location  of  the  maximum  heat 
transfer  rate  moves  away  from  the  step  as  the  free  stream  velocity 
increases. 

Abu-Mulaweh  et  al.  [40]  studied  experimentally  the  effect  of 
the  backward-facing  step  heights  on  the  turbulent  mixed  convec¬ 
tion  flow  along  a  3D  vertical  flat  plate.  The  experimental  setup,  step 
geometry,  boundary  conditions,  and  measurement  equipments 
were  identical  to  the  previous  work  done  by  the  same  authors  in 
2001  [39],  The  experiments  were  carried  out  for  step  heights  of 
0,  11,  and  22  mm,  while  the  free  stream  air  velocity,  u and  the 
temperature  difference  between  the  heated  walls  and  free  stream, 
AT,  were  fixed  at  0.41  m/s  and  30  °C,  respectively.  Measurements 
of  the  flow  field  were  carried  out  at  a  location  upstream  of  the 
step  (x  =  -5  cm)  at  the  midplane  (z= 0)  of  the  plate’s  width  for  pure 
forced,  pure  natural,  and  mixed  convection  flows.  It  was  found  that 
there  was  a  velocity  overshoots  in  the  mixed  convective  case  due 
to  the  buoyancy  assisting  near  to  the  heated  wall  upstream  of  the 
step  for  the  dimensionless  streamwise  velocity  distributions.  It  was 
observed  that  increasing  the  step  height  enhances  the  turbulence 
intensity,  which  causes  the  flow  to  become  turbulent  downstream 
of  the  step.  This  is  because  the  existence  of  a  backward-facing  step 
triggers  transition  from  the  re-laminarized  flow  upstream  of  the 
step  to  turbulent  flow  downstream  of  the  step.  Furthermore,  it  was 
found  that  as  the  step  height  increases,  the  turbulent  intensity  of 
both  streamwise  and  transverse  velocity  fluctuations,  the  intensity 
of  the  of  temperature  fluctuations,  and  heat  transfer  rate  down¬ 
stream  of  the  step  increase.  The  length  of  the  recirculation  region 
increases  with  increasing  the  step  height,  which  results  the  loca¬ 
tion  of  the  maximum  heat  transfer  rate  moves  away  from  the  step 
as  the  step  height  increases. 


4.  Fundamentals  of  nanofluids 

Heat  transfer  is  commonly  encountered  in  many  engineering 
systems  for  heating  and  cooling.  The  high  performance  cooling  is 
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essential  in  many  industrial  systems  such  as  nuclear  plants,  HVAC, 
transportation,  energy  production,  supply  to  electronics,  and  paper 
production.  It  is  also  important  in  consumer  products  such  as  com¬ 
puters,  car  engines,  power  electronics,  and  high-powered  lasers  or 
X-rays.  Traditional  fluids  are  used  as  coolant  and  their  heat  transfer 
increase  with  thermal  conductivity  increment.  The  development 
in  technologies  to  perform  desired  performance,  reliable,  and  min¬ 
imize  size  of  products  reveal  to  increase  in  heat  loads  which  in 
some  cases  it  exceeds  25  kW,  therefore  liquids  and  gases  have  lim¬ 
ited  thermal  performance  due  to  their  low  thermal  conductivity. 
Many  researchers  tried  to  increase  the  thermal  conductivity  of 
the  conventional  fluids  by  suspending  millimeter-  or  micrometer¬ 
sized  solid  particles  in  liquids  since  the  thermal  conductivity  of 
solids  are  higher  than  liquids.  Maxwell  [41]  made  the  first  effort 
using  millimeter-sized  particles  suspended  in  water,  but  due  to 
the  large  size  and  high  density  the  small  particles  settled  out  of  the 
suspension. 

Many  trails  have  been  done  to  disperse  solid  particles  in  liq¬ 
uids  since  Maxwell’s  time  and  all  of  them  observed  settlement  due 
to  the  large  size  and  density  of  the  particles.  In  1993,  Choi  con¬ 
ceived  the  novel  concept  of  dispersing  nanometer-sized  particles 
into  base  fluids  and  he  proposed  the  name  of  Nanofluids  in  1995 
[42].  Nanofluids  are  dilute  suspensions  of  nanoparticles  smaller 
than  lOOnm,  which  belong  to  a  new  type  of  functional  composite 
materials.  These  solid-liquid  composites  are  very  stable  and  show 
higher  thermal  conductivity  and  higher  convective  heat  transfer 
performance  than  traditional  liquids.  Nanofluids  have  small  size 
which  allows  them  to  have  less  particle  sedimentation,  erosion, 
passage  clogging,  and  pressure  drop  than  microfluids  or  fluids  with 
larger  scaled  suspensions. 

The  interest  of  studying  nanofluids  have  been  increased  since 
Choi’s  concept  is  presented,  and  many  researches  were  conducted 
in  different  fields  and  facilities  to  understand  the  structures,  func¬ 
tions,  mechanisms,  applications,  environmental  impacts,  and  other 
characteristics  of  nanofluids.  There  were  some  researchers  made 
review  papers  such  as  Yu  et  al.  [43],  Wang  and  Mujumdar  [44],  Li 
et  al.  [45],  Kakaq  and  Pramuanjaroenkij  [46],  Wen  et  al.  [47],  Das 
et  al.  [48],  Wong  and  Castillo  [49],  and  Nsofor  [50].  They  declared 
that  nanofluids  will  produce  a  new  generation  of  cooling  technol¬ 
ogy,  and  further  researches  is  needed  to  improve  the  stability  to 
avoid  cluster  phenomena  in  nanofluids. 

4.3.  Production  of  nanoparticles 

The  invention  and  improvement  of  modern  technologies  per¬ 
mits  to  fabricate  materials  with  nanometer  scale.  Many  types 
of  materials  have  been  used  to  make  nanoparticles  that  used  in 
nanofluids.  Nanoparticles  are  synthesized  by  two  different  pro¬ 
cesses  which  are  physical  process  and  chemical  process.  In  physical 
process,  mechanical  grinding  methods  using  grinding  balls,  and 
the  inert-gas-condensation  technique  were  utilized.  However, 
in  chemical  process  many  methods  are  used  such  as  chemical 
precipitation,  chemical  vapor  deposition,  micro-emulsions,  spray 
pyrolysis,  and  thermal  spraying.  Carbon  nano-tubes  (CNTs)  and 
fullerenes  are  commonly  formed  in  ordinary  flames  produced  by 
burning  methane,  ethylene,  and  benzene,  and  highly  irregular  in 
size  and  quality  in  nature.  These  types  of  nanoscaled  particles  are 
formed  by  different  processes  such  as  arc  discharge,  laser  ablation, 
chemical  vapour  deposition,  and  ball  milling. 

4.2.  Nanoparticle  material  types 

Certain  materials  have  been  used  to  produce  nanoparticles  and 
disperse  them  in  host  liquids  such  as  oxide  ceramics  (A1203,  CuO); 
nitride  ceramics  (AIN,  SiN);  carbide  ceramics  (SiC,  TiC);  metals 
(Ag,  Al,  Au,  Cu,  Fe);  semiconductors  (Si02,  Ti02);  single,  double, 


or  multi  wall  carbon  nanotubes  (SWCNT,  DWCNT,  MWCNT);  and 
composite  materials  such  as  nanoparticle  core-polymer  shell  com¬ 
posites.  Moreover,  new  materials  and  structures  are  desired  for  use 
in  nanofluids  where  the  solid-liquid  interface  is  doped  with  various 
molecules. 

4.3.  Host  liquid  types 

There  are  many  types  of  metallic  and  nonmetallic  liquids  have 
been  used  as  host  fluid  to  produce  nanofluids.  Metallic  liquids 
such  as  sodium  (644 1<)  have  higher  thermal  conductivity  than  non¬ 
metaflic  liquids  such  as  water  and  organic  liquids  (ethylene  glycol, 
oil).  Water  is  used  as  the  best  traditional  liquid  due  to  its  high 
thermal  conductivity,  abundance,  low  cost,  and  friendly  to  the  envi¬ 
ronment. 

5.  Preparation  of  nanofluids 

The  powder  form  nanoparticles  which  disperse  in  host  liq¬ 
uids  are  called  nanofluids.  Nanofluids  can  be  produced  by  two 
techniques;  the  two-step  (double-step)  method,  and  one-step 
(single-step)  method.  These  methods  have  been  utilized  using  dif¬ 
ferent  types  of  chemical  and  physical  techniques  to  make  sure  that 
the  solid-liquid  mixture  is  stable  to  avoid  agglomeration,  additional 
flow  resistance,  possible  erosion  and  clogging,  poor  thermal  con¬ 
ductivity,  and  poor  heat  transfer.  The  two-step  method  is  done 
by  producing  the  nanoparticles  powder  initially  as  introduced  in 
the  previous  section,  and  then  disperses  them  into  a  host  liquid. 
However,  in  one-step  method  the  nanoparticles  are  simultaneously 
made  and  directly  dispersed  into  the  base  fluid  [51  ].  It  is  noticed  in 
the  literature  that  nanofluids  with  oxide  nanoparticles  and  carbon 
nanotubes  are  produced  well  by  the  two-step  method,  while  it  is 
not  suitable  for  nanofluids  with  metallic  nanoparticles. 

5.3.  Techniques  of  nanofluids  production 

5.3.3.  Two-step  technique 

It  is  introduced  before  that  the  two-step  method  is  started 
by  nanoparticles  production  separately  and  afterward  dispersed 
into  a  base  fluid.  The  advantage  of  this  method  is  that  it  is  easily 
and  economically  produced.  On  the  other  hand,  the  disadvantage 
of  this  method  is  the  quick  agglomeration  of  individual  particles 
before  the  achievement  of  complete  dispersion  due  to  Vander 
Waals  attractive  forces  between  nanoparticles.  This  agglomera¬ 
tion  is  a  big  obstruction  to  achieve  high  heat  transfer  performance 
due  to  the  quick  settling  of  particles  out  of  the  base  fluids 
and  it  becomes  worse  as  the  volume  concentration  increases. 
The  agglomeration  is  not  only  a  problem  in  nanofluid  technol¬ 
ogy,  but  it  is  also  a  critical  issue  in  all  nanopowder  technology 
especially  during  drying,  storage,  and  nanoparticle  transportation 
stage. 

Many  different  methods  have  been  used  to  avoid  agglomeration 
using  two-step  process  to  initiate  the  move  towards  commercial¬ 
ization  by  facilitating  the  mass  production  of  nanofluids.  Two-step 
methods  such  as  stirrer,  ultrasonic  bath,  ultrasonic  disrupter,  and 
high  pressure  homogenizer  are  the  popular  methods  of  two-step 
technique.  Eastman  et  al.  [52],  Lee  at  al.  [53],  and  Wang  et  al.  [54] 
used  ultrasonic  equipment  to  disperse  A1203  in  the  base  fluid  with 
less  agglomeration.  Moreover,  Hong  et  al.  [55,56]  utilized  ultrasonic 
equipment  produced  Fe  nanofluids  by  mixing  Fe  nanocrystalline 
powder  in  ethylene  glycol.  Ultrasonic  equipment  was  also  used 
to  produce  Ti02 -water  nanofluid  as  introduced  in  Murshed  et  al. 
[57].  Metallic  nanofluids  were  produced  by  using  different  two-step 
methods,  but  they  were  not  successful  due  to  the  high  agglom¬ 
eration.  There  are  some  other  techniques  that  are  used  to  attain 
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the  stability  of  the  nanoparticles  in  the  base  fluids  using  the  two- 
step  methods  such  as  controlling  pH,  or  adding  surfactants  and 
dispersant.  These  techniques  are  explained  in  detail  in  the  next 
section. 

5.1.2.  One-step  technique 

One-step  technique  is  the  process  of  producing  nanofluids 
where  the  nanoparticles  are  simultaneously  made  and  directly  dis¬ 
persed  into  the  base  fluid.  This  method  is  preferable  to  produce 
nanofluids  containing  high  thermal  conductivity  metals  to  avoid 
erosion  and  oxidation  of  particles.  The  advantage  of  this  process 
lies  on  minimizing  of  nanoparticles  agglomeration.  This  behavior 
increases  the  stability  of  the  suspensions  and  uniform  dispersion 
in  the  host  liquids.  The  disadvantage  of  one-step  technique  is  the 
limit  of  quantity  of  the  production  due  to  the  slow  of  the  production 
process,  low  concentration  of  nanoparticles,  and  the  high  cost. 

Different  methods  have  been  used  to  reduce  the  time  and  cost 
using  evaporation,  physical  or  chemical  one-step  methods.  The 
first  approach  of  one-step  evaporation  method  is  direct  evapora¬ 
tion  method  which  called  Vacuum  Evaporation  onto  a  Running  Oil 
Substrate  (VEROS)  technique  developed  by  Akoh  et  al.  [58],  This 
method  was  created  to  produce  nanoparticles,  but  it  was  found 
that  it  is  difficult  to  separate  the  particles  from  fluids  to  produce 
dry  nanoparticles.  Later,  Wagener  et  al.  [59]  proposed  a  modified 
VEROS  technique  using  high  pressure  magnetron  sputtering  to  pro¬ 
duce  metallic  nanofluids  such  as  Ag  and  Fe.  Furthermore,  Eastman 
et  al.  [52]  modified  VEROS  technique  by  using  direct  condensed 
Cu  vapour  to  form  nanoparticles  and  in  contact  with  flowing  with 
low-vapor  pressure  ethyl  glycol  (EG).  The  well  dispersed  of  Cu 
nanoparticles  in  ethyl  glycol  enhanced  the  thermal  conductivity  of 
the  base  fluid  by  up  to  40%  at  the  particle  volume  concentration  of 
0.3  vol.%.  For  one-step  physical  method;  a  submerged  arc  nanopar¬ 
ticles  synthesis  is  used  in  [60-62]  to  produce  nanofluids  of  Cu,  CuO, 
and  TiC>2.  In  this  method,  nanoparticles  are  produced  by  heating 
the  solid  material  from  an  electrode  in  the  form  of  arc  sparkling 
and  condensed  into  liquid  in  a  vacuum  chamber  to  form  nanofluid. 
Finally,  the  one-step  chemical  method  is  reported  by  Zhu  et  al.  [63]. 
They  produced  Cu  nanoparticles  of  less  than  20  nm  and  dispersed 
in  (EG)  by  the  reduction  of  a  copper  salt  by  sodium  hypophosphite 
under  microwave  irradiation.  A  protective  polymer  and  stabilizer 
called  polyvinylpyrrolidone  was  added  so  that  it  inhibited  agglom¬ 
eration.  Wei  et  al.  [64]  developed  a  chemical  solution  method  (CSM) 
to  synthesize  Cu20  nanofluids  by  suspensions  of  cuprous-oxide 
(Cu20)  in  water.  Their  experimental  studies  showed  a  substantial 
thermal  conductivity  enhancement  of  up  to  24%  with  their  synthe¬ 
sized  nanofluids.  Among  all  the  three  methods;  chemical  process 
is  the  fastest  and  can  be  produced  in  large  quantities  but  it  is  still 
not  in  commercial  quantity  as  in  two-step  process. 

5.2.  Stability 

The  stability  of  the  suspensions  in  nanofluids  is  the  most  criti¬ 
cal  issue  for  enhancing  heat  transfer.  This  is  because  nanoparticles 
tend  to  aggregate  with  the  time  elapsed  for  its  high  surface-activity, 
which  results  not  only  the  settlement  and  clogging  of  microchan¬ 
nels  but  also  the  decreasing  of  thermal  conductivity  of  nanofluids. 
Nanofluids  are  considered  to  be  stable  when  the  concentration 
or  particle  size  of  supernatant  particles  keeps  constant.  There  are 
some  methods  have  been  used  to  determine  the  stability  of  nanoflu¬ 
ids  such  as  sedimentation  method  [65,66],  Zeta  potential  method 
and  absorbency  method  [67].  Sedimentation  method  is  the  most 
simple  and  reliable  method.  It  is  basically  obtained  by  special  appa¬ 
ratus  which  functioning  to  measure  the  variation  of  concentration 
or  particle  size  of  supernatant  particle  with  sediment  time.  More¬ 
over,  another  method  using  an  optical  sedimentation  photograph 
to  monitor  the  stability  of  nanofluids  was  used  but  its  disadvantage 


is  the  long  time  of  the  observation  compared  to  other  methods.  On 
the  other  hand,  Zeta  potential  method  has  a  limitation  of  the  viscos¬ 
ity  and  concentration  of  nanofluids.  A  UV-vis  spectrophotometer 
was  used  to  detect  the  stability  of  nanofluids.  The  variation  of  super¬ 
natant  particle  concentration  of  nanofluids  with  sediment  time 
can  be  obtained  by  the  measurement  of  absorption  of  nanofluids 
because  there  is  a  linear  relation  between  the  supernatant  nanopar¬ 
ticle  concentration  and  the  absorbance  of  suspended  particles.  It  is 
declared  that  the  most  important  factors  affecting  the  stability  of 
suspensions  were  the  nanoparticles  concentration,  additives,  vis¬ 
cosity  of  base  liquid  and  pH  value.  The  variety,  diameter,  density  of 
nanoparticle  and  ultrasonic  vibration  also  influence  the  stability  of 
nanofluids.  In  addition,  the  particle  morphology,  chemical  structure 
of  particles  and  base  fluids  affect  the  stability  as  well.  It  is  observed 
that  the  stability  of  nanofluids  do  not  exceed  more  than  few  weeks 
or  months  which  need  to  have  new  techniques  to  improve  the  time 
of  nanofluids  stability. 

6.  Experimental  investigations  on  nanofluids 

6.1.  Thermal  conductivity 

It  is  reported  that  many  experimental  studies  have  been  done  on 
thermal  conductivity  of  nanofluids  because  it  is  a  primary  assess¬ 
ment  of  the  heat  transfer  performance  of  nanofluids.  The  physical 
mechanism  for  thermal  conductivity  enhancement  of  nanofluids  is 
not  well  understood  up  to  date.  Maxwell  [68]  was  one  of  the  first 
investigators  who  suspended  particles  thermal  conductivity.  He 
used  a  very  dilute  suspension  of  spherical  particles  by  ignoring  the 
interactions  between  the  particles,  and  made  an  equation  that  can 
be  applied  only  to  mixtures  with  low  particle  volume  concentra¬ 
tions.  The  gap  of  particles  interactions  attracts  many  researchers  to 
extend  or  improve  the  Maxwell  equation  with  various  factors  that 
affect  thermal  conductivity  such  as  particle  shape  [69-77],  parti¬ 
cle  dispersion  and  distribution  [78],  volume  concentration  [79-82], 
particle-shell  structure  [69,83-88],  and  contact  resistance  [89,90]. 
Even  though  these  studies  tried  to  improve  equations  to  predict 
thermal  conductivity  but  the  comparison  between  the  experimen¬ 
tal  data  and  theoretical  predictions  for  nanofluids  is  generally  not 
satisfactory.  New  techniques  of  mechanism  have  been  used  to 
improve  the  prediction  of  nanoscale  such  as  nanoparticle-matrix 
interfacial  layer  [91-94],  nanoparticle  Brownian  motion  [95-99], 
and  nanoparticle  cluster/aggregate  [100,101]. 

6. 1. 1.  Thermal  conductivity  measurement  methods 

There  are  many  methods  used  to  measure  nanofluids  thermal 
conductivity  such  as  transient  hot-wire  (THW),  temperature  oscil¬ 
lation  (TO)  and  steady-state  (SS)  methods  [55-57,61,102-105]. 
Transient  hot-wire  method  is  well  established  as  the  most  accu¬ 
rate,  reliable  and  robust  measurement  technique  for  the  thermal 
conductivity  of  nanofluids,  and  so  it  is  utilized  in  most  of  thermal 
conductivity  measurements.  The  transient  hot-wire  is  functioned 
as  both  heater  and  thermometer  and  many  modifications  have  been 
proposed  such  as  coating  the  hot  wire.  In  oscillation  method;  the 
electrical  components  are  removed  from  the  test  sample  and  per¬ 
form  as  purely  thermal  as  reported  by  [106-108],  so  that  the  ion 
movement  is  not  affecting  the  measurements.  The  steady-state 
method  was  used  by  [54,109]  to  measure  the  effective  thermal 
conductivity  of  A1203  and  CuO  nanofluids  with  an  improvement 
of  12%  in  thermal  conductivity  enhancement  at  3%  volume  fraction 
of  nanoparticles. 

6. 1.2.  Effects  of  thermal  conductivity  parameters 

Many  experimental  studies  have  been  done  by  different 
researchers  for  examining  and  improving  the  thermal  conduc¬ 
tivity  enhancement  of  certain  types  of  nanofluids  using  different 
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methods  and  parameters.  They  calculated  the  thermal  conductivity 
enhancement  by  using  the  ratio  of  the  thermal  conductivity  of  the 
nanofluid  over  its  base  fluid  thermal  conductivity.  It  is  noticed  that 
most  of  the  studies  did  not  concerned  or  noticed  the  particle  size 
distribution.  In  addition,  the  effect  of  the  particles  agglomeration 
in  the  base  fluid  is  also  unknown. 

It  is  declared  that  the  thermal  conductivity  of  nanofluid  is  sen¬ 
sitive  to  certain  parameters  which  are  responsible  for  constant 
nanofluids  such  as  particle  volume  concentration,  particle  size, 
particle  material,  particle  shape,  base  fluid  material,  temperature, 
additives  and  surfactants,  and  acidity.  These  parameters  will  be 
introduced  in  the  following  section  separately  base  on  the  previous 
experimental  investigations  and  findings. 

6.1. 2.1.  Effects  of  particle  volume  concentration.  One  of  the  most 
important  parameters  is  volume  concentration  or  in  other  words 
volume  fraction.  It  is  the  ratio  of  the  volume  of  nanoparticles  to 
the  base  fluid  volume.  As  declared  previously  that  solids  have 
higher  thermal  conductivity  than  liquids,  so  as  the  percentage 
of  nanoparticles  increases,  the  thermal  conductivity  increases  as 
well.  Unfortunately,  this  is  not  working  for  all  cases;  it  is  declared 
that  many  colloidal  and  biological  suspensions  show  strong  non- 
Newtonian  behavior  but  nanofluids  with  volume  concentration  not 
more  than  5%  exhibit  Newtonian  behavior.  At  higher  particle  vol¬ 
ume  concentrations,  the  increase  in  enhancement  is  expected  to 
diminish  or  even  reverse,  and  nanofluid  viscosity  would  increase. 
Three  methods  have  been  used  by  Wang  et  al.  [54]  to  measure  the 
viscosity  of  nanofluids  but  they  did  not  observe  any  non-Newtonian 
effects.  They  used  Al2C>3-water  nanofluid  and  found  a  30%  increase 
in  viscosity  when  compared  with  pure  water  at  3%  volume  con¬ 
centration  of  the  particles.  However,  Pak  and  Cho  [110]  measured 
the  viscosity  to  be  much  higher.  There  are  many  researches  done 
on  particle  volume  concentration  [42,54, 1 02, 1 06, 1 09, 1 1 1  - 1 26]  for 
different  particle  size,  particle  shape,  particle  material,  base  fluid 
material,  additive,  and  acidity  level.  They  reported  that  at  low  vol¬ 
ume  concentration  (<5%)  the  thermal  conductivity  ratios  increase 
linearly  with  volume  concentration.  Furthermore,  it  is  noticed  that 
metallic  nanofluids  and  either  single-  or  multi-wall  nanotubes 
enhance  thermal  conductivity  with  low  volume  concentration 
compared  to  oxides  and  nonmetallic  nanofluids. 

6.1. 2.2.  Effects  of  particle  size.  The  size  of  particle  plays  a  significant 
role  in  thermal  conductivity  and  heat  transfer  enhancement  in  base 
fluids.  As  introduced  before  that  different  suspension  sizes  of  mil¬ 
limeter  and  micrometer  particles  in  host  liquid  have  been  used  to 
increase  the  thermal  conductivity,  but  the  particles  agglomerate 
quickly  and  settled  out  of  the  liquid.  Nanoparticles  are  formed  to 
decrease  the  size  of  the  particles  so  that  the  time  of  sedimenta¬ 
tion  decreases  to  reach  in  some  cases  to  more  than  few  weeks  or 
months.  The  size  of  nanoparticles  are  in  the  range  of  l-100nm  and 
as  the  size  decreases  the  thermal  conductivity  increases  due  to  the 
increase  of  the  relative  surface  area  [127],  The  much  larger  relative 
surface  area  of  nanoparticles,  should  not  only  significantly  improve 
thermal  conductivity  capabilities,  but  also  should  increase  the  sta¬ 
bility  of  the  suspensions.  In  addition,  it  is  reported  that  as  the  size  of 
nanoparticle  reduced,  the  Brownian  motion  will  be  induced  [128]. 
Furthermore,  there  are  some  researches  indicate  that  the  thermal 
conductivity  enhanced  with  nanoparticles  increment  to  Ri60nm 
and  vice  versa  for  greater  values  [53,114].  The  size  of  nanopar¬ 
ticles  in  the  host  liquid  is  too  important  to  the  modern  science 
and  technology.  It  can  be  applied  in  small  size  systems  such  as 
micro-channels  and  nanosystems  where  the  clog  and  erosion  is 
prohibited. 

6.1.23.  Effects  of  particle  material.  There  are  many  types  of  materi¬ 
als  have  been  utilized  to  produce  nanofluids  such  as  oxide  ceramics 


(AI2O3,  CuO);  nitride  ceramics  (AIN,  SiN);  carbide  ceramics  (SiC, 
TiC);  metals  (Ag,  Al,  Au,  Cu,  Fe);  semiconductors  (SiC>2,  Ti02);  sin¬ 
gle,  double,  or  multi  wall  carbon  nanotubes  (SWCNT,  DWCNT, 
MWCNT);  and  composite  materials.  It  was  found  that  metallic 
and  allotrops  of  carbon  nanoparticles  (such  as  diamond)  increase 
the  thermal  conductivity  of  a  host  fluid  more  than  other  types  of 
nanoparticles  having  the  same  volume  concentration. 

6.1. 2.4.  Effects  of  particle  shape.  Thermal  conductivity  enhance¬ 
ment  in  nanofluids  is  affected  by  the  geometrical  shape  of  the 
particles.  Different  shapes  such  as  spherical  and  cylindrical  (fiber, 
rod,  and  tube)  nanoparticles  have  been  used  in  different  host  liquids 
to  identify  the  enhancement  of  nanofluids  thermal  conductivity.  It 
is  reported  that  cylindrical  nanoparticles  are  better  than  spherical 
shape  in  enhancing  the  thermal  conductivity  of  nanofluids  [57,129]. 
This  is  due  to  the  elongated  particles  and  the  large  aspect  ratio  of 
the  surface  area  to  volume  that  conducts  heat  through  the  fluid. 

6.1. 2.5.  Effects  of  base  fluid  material.  Many  host  liquids  have  been 
used  to  produce  nanofluids  such  as  aqueous  and  organic  liquids 
(ethylene  glycol,  and  oils)  to  enhance  thermal  conductivity.  It  is 
found  that  the  thermal  conductivity  enhancement  increases  as  the 
thermal  conductivity  of  the  base  fluid  decreases  [114,115,130]. 

6.1. 2.6.  Effects  of  temperature.  The  nature  of  the  nanofluid  ther¬ 
mal  conductivity  is  sensitive  to  temperature.  Many  studies  have 
been  done  for  the  effect  of  temperature  on  thermal  conductiv¬ 
ity.  They  reported  that  the  thermal  conductivity  enhanced  as  the 
temperature  increases  due  to  the  motion  of  the  nanoparticles 
[106,109,127,131].  In  opposite,  Masuda  et  al.  [Ill]  found  that  the 
thermal  conductivity  enhancement  decreases  as  the  temperature 
increases. 

6.1. 2. 7.  Effects  of  additives  and  surfactants.  It  is  declared  that  certain 
types  of  additive  have  been  used  in  nanofluids  to  keep  nanoparti¬ 
cles  in  suspension  form  and  prevent  agglomeration  for  improving 
the  thermal  conductivity  enhancement.  The  additives  types  and 
concentration  is  added  depend  on  the  type  of  both  nanoparticles 
and  base  fluids.  Eastman  et  al.  [113]  used  thloglycolic  acid  in  Cu- 
ethyl  glycol,  Assael  et  al.  [132]  used  CTAB  and  Nanosperse  AQ  for 
MWCNT-water,  Ding  et  al.  [122]  gum  Arabic  for  MWCNT-water, 
and  other  types  of  acidic  fluids  such  as  Oleic  acids,  etc. 

6. 1.2.8.  Effects  of  Acidity  (PH).  Limited  studies  have  been  published 
on  the  effect  of  fluid  acidity  on  the  thermal  conductivity  enhance¬ 
ment  of  nanofluids  [117,126].  But  the  general  trend  is  that  acidity 
increases  the  thermal  conductivity  enhancement. 

6.1.3.  Mechanisms  of  thermal  transport  in  nanoscale 
6.13.1.  Brownian  motion.  Brownian  motion  is  one  of  the  param¬ 
eters  that  affect  the  flow  and  heat  transfer  of  nanofluids.  It  is 
basically  the  random  dynamic  mode  of  particles  in  a  liquid  where 
the  particles  collide  between  each  other.  It  was  determined  that 
the  thermal  conductivity  increases  with  an  increase  in  tempera¬ 
ture,  but  also  it  was  shown  that  nanofluids  compose  of  smaller 
particles  experienced  greater  enhancement  than  with  larger  par¬ 
ticles.  Since  temperature  represents  the  overall  kinetic  energy  of 
the  particles,  an  increase  in  temperature  will  cause  increase  in 
the  particles  motion.  It  is  easier  for  smaller  particles  to  move; 
therefore,  smaller  particles  will  display  a  higher  level  of  Brown¬ 
ian  motion  than  larger  particles.  It  is  outlined  that  the  movement 
of  nanoparticles  due  to  Brownian  motion  was  too  slow  in  trans¬ 
porting  heat  through  a  fluid  and  it  can  be  ignored  [133-136],  To 
travel  from  one  point  to  another,  a  particle  moves  a  large  dis¬ 
tance  over  many  different  paths  in  order  to  reach  a  destination 
that  may  be  a  short  distance  from  the  starting  point.  Therefore, 
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the  random  motion  of  particles,  no  matter  how  agitated  or  ener¬ 
getic  they  may  be  not  a  key  factor  in  the  improvement  of  heat 
transfer.  On  the  other  hand,  there  are  some  studies  developed  a 
dynamic  model  that  takes  into  account  convection  heat  transfer 
induced  by  Brownian  nanoparticles.  They  reported  that  Brownian 
motion  of  nanoparticles  at  the  molecular  and  nanoscale  levels  is  a 
main  mechanism  controlling  the  thermal  conductivity  of  nanoflu¬ 
ids  [96,98,137], 

6.13.2.  Liquid  layering  on  the  nanoparticle-liquid  interface.  There 
are  some  studies  declared  that  there  is  a  liquid  layering  on  the 
nanoparticles,  which  improves  the  enhancement  of  nanofluid  heat 
transfer  properties.  The  thickness  and  thermal  conductivity  of  the 
nanolayer  are  not  known  yet,  but  the  liquid  molecules  close  to  a 
solid  surface  have  been  proven  by  Yu  et  al.  [138]  to  form  layers. 
Ren  et  al.  [94]  indicated  that  the  thermal  conductivity  increases 
with  liquid  layer  increment.  They  assumed  that  the  thermal  con¬ 
ductivity  of  the  layer  would  be  somewhere  between  the  thermal 
conductivity  of  the  bulk  fluid  and  the  nanoparticle.  On  the  other 
hand,  Keblinski  [135]  assumed  that  the  thermal  conductivity  of  the 
interfacial  liquid  is  the  same  as  that  of  the  solid  and  showed  that 
the  resultant  large  effective  volume  of  the  particle-layered-liquid 
structure  would  enhance  the  thermal  conductivity.  Ren  et  al.  [94] 
also  found  that  as  the  nanoparticles  increased  in  size,  the  effects  of 
the  liquid  layering  became  weaker.  Yu  and  Choi  [92]  determined 
that  the  solid-liquid  interfacial  layers  in  the  nanofluid  plays  an 
important  role  in  the  enhanced  thermal  conductivity  especially 
when  the  particle  diameter  is  less  than  10  nm. 

6.1.3.3.  Nature  of  heat  transport  in  nanoparticles.  Heat  is  transferred 
via  conduction  implies  that  heat  is  transferred  internally,  by  vibra¬ 
tions  of  atoms  and  molecules.  The  vibrations  of  the  atoms  which  are 
connected  together  as  give  rise  to  the  vibrations  of  the  whole  crys¬ 
tal,  that  is,  to  lattice  vibrations.  The  energy  of  the  whole  vibrating 
system  is  quantized,  and  the  quantum  of  thermal  energy  absorbed 
or  emitted  by  an  atom  is  called  a  phonon.  Therefore,  a  phonon  is 
essentially  a  quantized  mode  of  vibration  occurring  in  rigid  crys¬ 
tal  lattice  and  plays  a  major  role  in  material’s  thermal  conduction. 
Keblinski  [135]  demonstrated  that  the  phonon  mean  free  bath  is 
much  shorter  in  the  liquid  than  in  the  particle  and  would  be  effec¬ 
tive  as  the  distance  between  the  particles  and  the  thickness  of  the 
layered  liquid  are  small.  A  greater  phonon  density  exists  in  hot 
regions  of  a  crystal  than  in  cooler  regions.  Therefore,  heat  conduc¬ 
tion  is  essential  due  to  the  diffusion  of  phonons,  and  because  the 
temperature  gradient  changes  from  hot  to  cold  regions.  Electrons 
can  also  carry  heat;  thus,  metals  have  many  free  electrons  which 
move  around  randomly,  heat  can  be  transferred  from  one  part  of  the 
metal  to  another  quite  effectively.  That  is  why  metals  are  generally 
very  good  conductors  of  heat. 

6. 3.3.4.  Effects  of  nanoparticles  clustering.  Clustering  of  particles  is 
one  of  the  major  problems  that  affects  the  enhancement  of  ther¬ 
mal  conductivity  by  creating  path  of  lower  thermal  resistance. 
It  is  noticed  that  as  the  particles  agglomerate  increases,  the  size 
increases  and  becomes  denser  where  they  clustered  down,  and 
end  with  low  thermal  conductivity  enhancement.  To  avoid  or  min¬ 
imize  the  particles  clustering,  nanofluids  should  consist  of  low 
volume  concentration  and  small  particles  size.  The  effective  vol¬ 
ume  of  a  cluster,  which  is  the  volume  from  which  other  clusters 
are  excluded,  can  be  much  larger  than  the  physical  volume  of  the 
particles  [135,136].  The  effective  volumetric  concentration  of  the 
particles,  inside  a  cluster,  is  larger  than  the  volume  of  the  solid 
phase,  since  within  such  cluster  heat  can  move  very  rapidly. 

6. 3.3.5.  Thermophoresis.  Thermophoresis  refers  to  the  motion  of 
colloidal  particles  in  response  to  a  temperature  gradient  which 


can  be  clarified  by  applying  the  kinetic  theory.  The  high  energy 
molecules  in  the  hot  region  of  the  liquid  impinge  on  the  particles 
with  greater  momentum  than  do  molecules  coming  from  the  cold 
region,  thus  leading  to  the  migration  of  the  particles  in  the  direction 
opposite  to  the  temperature  gradient.  This  mechanism  has  been 
found  to  be  negligible,  primarily  because  of  the  very  small  particle 
migration  velocity  [91,136]. 

6. 3.3.6.  Reduction  in  thermal  boundary  layer  thickness.  A  few  num¬ 
bers  of  researches  have  been  done  for  studying  the  thermal 
boundary  layer  such  as  in  Ding  et  al.  [139].  They  mentioned 
that  a  reduction  in  the  thermal  boundary  layer  thickness  may 
be  a  mechanism  that  causes  heat  transfer  enhancements  in 
nanofluids.  It  is  still  needed  to  have  more  research  on  this  mech¬ 
anism  to  understand  its  effects  in  nanofluids  convective  heat 
transfer. 

6.2.  Viscosity 

Nanofluids  are  expected  to  be  used  under  flow  conditions  and 
the  flow  of  suspension  which  are  sometimes  drastically  different 
from  that  of  most  common  heat  transfer  fluids  that  have  Newtonian 
characteristics,  so  it  is  essential  to  have  the  rheological  properties  of 
nanofluid  to  use  practically.  Rheological  properties  can  provide  the 
knowledge  on  the  microstructure  under  both  static  and  dynamic 
conditions,  which  are  important  to  understand  the  mechanism  of 
heat  transfer  enhancement  using  nanofluids.  Fluids  have  frictional 
forces  between  the  molecules  and,  therefore,  they  display  a  certain 
flow  resistance  which  can  be  measured  as  viscosity.  The  fluids  are 
either  Newtonian  or  non-Newtonian,  depends  on  their  viscosity.  It 
has  been  demonstrated  by  Newton  that  the  shear  force  acting  on 
a  liquid  is  proportional  to  the  resulting  flow  velocity.  If  the  viscos¬ 
ity  of  the  fluid  remains  constant  with  an  increase  in  shear  rate,  it  is 
called  Newtonian  fluid,  but  if  the  viscosity  changes  with  an  increase 
in  shear  rate,  it  is  called  non-Newtonian  fluid.  The  shear  viscosity 
of  samples  displaying  shear-thickening  behavior  is  dependent  on 
the  degree  of  shear  load.  However,  the  viscosity  increases  with  an 
increase  in  shear  stress.  It  is  declared  that  the  viscosity  of  nanoflu¬ 
ids  is  mainly  dependent  on  volume  concentration,  nanoparticle 
size,  and  temperature  [140-148].  Nanofluids  with  high  concen¬ 
tration  and  large  size  nanoparticles  have  high  particle  interaction 
and  may  result  in  particles  agglomeration  and  thus  increase  the 
flow  resistance  and  decrease  the  heat  transfer  enhancement.  Par¬ 
ticle  shape  plays  an  important  role  since  during  the  shear  process, 
particles  rotate  as  they  move.  Cube-shape  particles  take  up  more 
volume  when  rotating  than  spherical  particles  and  hence  less  free 
volume  is  available  for  the  liquid  between  the  particles.  In  addi¬ 
tion,  it  was  reported  that  the  viscosity  decreases  with  temperature 
increment. 

6.3.  Heat  capacity 

The  heat  capacity  of  a  substance  is  the  amount  of  heat  required 
to  change  its  temperature  by  one  Kelvin,  and  has  units  of  joule  per 
Kelvin  (J/K)  in  the  SI  system.  The  specific  heat  capacity  is  essen¬ 
tially  the  heat  capacity  per  unit  mass,  usually  per  gram  of  material 
(J/gK).  The  specific  heat  capacities  of  nanofluids  are  different  from 
that  of  base  fluid  and  increase  with  the  size  and  volume  concen¬ 
tration  of  nanoparticles  decrement.  The  high  specific  interfacial 
area  of  nanoparticle  can  absorb  liquid  molecules  to  its  surface  and 
form  liquid  layers,  which  will  reversely  constrain  nanoparticle  and 
turns  its  free  boundary  surface  atoms  to  be  non-free  interior  atoms. 
This  effect  will  also  increase  the  specific  heat  capacity  of  nanofluid 
[149-153], 
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7.  Applications  of  nanofluids 

Nanofluids  could  be  used  as  a  super  coolant  in  many  industrial 
and  consumer  products  such  as  nuclear  reactors,  car  engines,  radi¬ 
ators,  computers,  X-rays,  and  many  other  applications.  It  is  called 
super-coolant  because  it  can  absorb  heat  more  than  any  traditional 
fluids,  so  it  can  minimize  the  size  of  system  and  increase  its  per¬ 
formance.  Most  of  the  researches  were  concerned  on  studying  the 
properties,  parameters,  and  heat  transfer  enhancement,  but  not 
much  research  done  for  commercial  systems.  Some  researchers 
tried  to  apply  nanofluids  in  real  systems  and  observe  their  flow 
and  heat  transfer.  Tzeng  et  al.  [154]  studied  the  effect  of  nanofluids 
when  used  as  engine  coolants.  CuO  (4.4wt%)  and  AI2O3  (4.4wt%) 
nanoparticles  and  antifoam  were  individually  mixed  with  auto¬ 
matic  transmission  oil,  and  experiment  them  by  a  real  drive  of 
four-wheel-drive  (4WD)  transmission  system.  The  experimental 
results  showed  that,  the  antifoam-oil  provided  the  highest  tem¬ 
perature  distribution  in  rotary  blade  coupling  and,  accordingly,  the 
worst  heat  transfer  effect,  and  CuO-oil  provided  the  lowest  tem¬ 
perature  distribution  both  at  high  and  low  rotating  speed  and, 
accordingly,  the  best  heat  transfer  effect.  Chein  and  Huang  [155] 
used  Cu  nanofluid  with  various  concentrations  to  study  silicon 
microchannel  heat  sink  performance.  They  found  that  nanofluid 
cools  better  than  pure  water  and  observed  no  extra  pressure  drop 
because  of  small  particle  size  and  low  particle  volume  fraction. 
Koo  and  Kleinstreuer  [156]  used  two  types  of  nanofluids  (i.e.,  CuO 
nanospheres  at  low  volume  concentrations  in  water  and  in  ethy¬ 
lene  glycol)  for  the  conjugate  heat  transfer  problem  for  microheat 
sinks.  The  effect  of  Brownian  motion  on  the  effective  fluid  viscos¬ 
ity  was  considered  and  found  to  be  less  significant  than  that  on 
the  effective  thermal  conductivity.  Many  other  industrial  and  con¬ 
sumer  companies  have  their  own  researches  using  nanofluids  such 
as  medical  equipments,  computer  manufacturing,  and  gas  and  oil 
[157], 

8.  Safety 

Up  to  date  the  behavior,  parameters  and  mechanisms  of 
nanofluids  are  not  clearly  understood,  and  because  of  these  vari¬ 
ables  scientists  still  cannot  accurately  predict  how  nano-materials 
will  affect  living  organisms.  Now,  it  is  just  clear  that  the  biolog¬ 
ical  activity  and  biokinetics  of  nanoparticles  are  different  from 
larger  particles,  and  that  they  depend  on  many  parameters.  These 
parameters  can  modify  cellular  uptake,  protein  binding,  transloca¬ 
tion  from  portal  of  entry  to  target  site,  and  the  possibility  of  tissue 
injury.  Several  studies  have  demonstrated  that  nanoparticles  are 
more  toxic  than  micro-particles  of  the  same  material. 

9.  Environmental  effects 

Nanofluids  are  observed  to  be  an  environmental  friend  due  to  its 
high  heat  transfer  enhancement.  The  usage  of  nanofluids  decreases 
the  consumption  of  fuel  in  engines,  and  electricity  which  are  the 
highest  pollutants  that  threat  our  earth  and  increase  its  temper¬ 
ature.  It  is  noticed  that  if  nanofluids  improve  chiller  efficiency  by 
1%;  it  will  save  320  billion  kwh  of  electricity  or  equivalent  to  5.5 
million  barrels  of  oil  per  year  just  in  USA  alone. 

10.  Convective  heat  transfer  of  nanofluids 

Convective  heat  transfer  is  the  phenomena  of  the  macroscopic 
motion  of  the  fluid  relative  to  the  surface  due  to  the  heat  trans¬ 
fer  between  them.  The  surface  could  be  a  solid  wall  or  an  interface 
with  another  liquid.  Convective  heat  transfer  is  divided  into  two 
major  types,  which  are  natural  convective  heat  transfer  where  fluid 


motion  is  induced  by  buoyancy,  and  forced  convective  heat  transfer 
where  fluid  is  forced  to  flow  through  a  boundaiy  region.  Convective 
heat  transfer  of  a  medium  depends  on  its  thermophysical  proper¬ 
ties  such  as,  thermal  conductivity,  specific  heat  capacity,  viscosity, 
density,  and  thermal  expansion  coefficient. 

10.  J.  Natural  convection 

The  natural  convection  studies  using  nanofluids  are  very  few. 
Even  though  the  research  is  limited,  however  the  available  results 
are  diverged.  Putra  et  al.  [158],  and  Wen  and  Ding  [120]  made 
an  experimental  investigation  for  natural  convection  using  AI2O3 
and  TiC>2  nanofluids,  respectively.  They  found  that  the  presence  of 
nanoparticles  in  water  systematically  decreased  the  natural  con¬ 
vective  heat  transfer  coefficient.  Furthermore,  they  declared  that 
the  decrease  was  due  to  the  effects  of  particle/fluid  slip  and  sedi¬ 
mentation  of  nanoparticles.  In  contrast,  Khanafer  et  al.  [159]  used  a 
numerical  technique  and  predicted  that  nanofluids  enhanced  nat¬ 
ural  convective  heat  transfer.  It  is  needed  to  do  more  research  on 
natural  convection  to  get  a  clear  view  and  well  understanding  for 
such  phenomena. 

1 0.2.  Forced  convection 

Forced  convection  studies  using  nanofluids  are  also  limited  and 
they  have  diverged  results  from  each  other.  Majority  of  forced  con¬ 
vective  heat  transfer  using  nanofluids  found  that  there  is  convective 
heat  enhancement  such  as  Li  and  Xuan  [160],  Xuan  and  Li  [161], 
Lee  and  Choi  [162],  Jang  and  Choi  [163],  and  Xuan  and  Roetzel 
[164].  There  are  some  studies  reported  that  heat  transfer  can  be 
enhanced  under  certain  conditions  (such  as  low  volume  concen¬ 
tration  for  decreasing  the  viscosity)  by  Pak  and  Cho  [165],  Lee  and 
Mudawar  [166],  and  Chein  and  Chuang  [167].  A  very  few  studies 
show  a  decrease  in  the  convective  heat  transfer  coefficient  when 
nanoparticles  are  added  to  the  base  liquids  such  as  Yang  et  al.  [  1 68  ]. 

10.3.  Convection  of  nanofluids  in  BFS 

It  is  noticed  in  the  literature  that  there  are  few  experimen¬ 
tal  and  numerical  studies  have  been  done  to  investigate  the  flow 
and  heat  transfer  of  nanofluids  in  different  geometries  with  dif¬ 
ferent  regimes.  Engineering  geometries  such  as  tubes  and  annuli 
were  the  first  geometries  that  have  been  used  for  investigating 
nanofluids  enhancement.  Most  of  the  studies  show  that  nanoflu¬ 
ids  have  higher  thermal  transport  than  any  ordinary  fluids  and 
their  enhancement  increase  or  decrease  depends  on  the  change 
of  the  parameters  and  characteristics  which  introduced  in  the  pre¬ 
vious  sections.  Other  studies  have  been  done  to  apply  nanofluids 
for  microchannels,  different  geometrical  shapes,  and  systems.  The 
results  show  an  improvement  in  thermal  transport  enhancement. 

To  the  best  knowledge  of  the  authors,  there  is  only  one  paper 
investigated  the  flow  and  heat  transfer  influence  of  nanofluids 
through  a  duct  with  backward  facing  step.  Abu-Nada  [169]  stud¬ 
ied  the  flow  and  heat  transfer  influences  over  a  backward  facing 
step  (BFS)  using  nanofluids  numerically.  Different  volume  fractions 
and  types  of  nanoparticles  were  used  in  the  base  fluid  (water)  and 
flowed  through  BFS  channel  with  an  expansion  ratio  ER  =  2.  The 
Reynolds  number  and  nanoparticles  volume  fraction  used  were 
in  the  range  of  200  <  Re  <600  and  0<ip<0.2,  respectively.  The 
Prandtl  number  of  the  base  fluid  was  kept  constant  at  6.2  for  the 
five  types  of  nanoparticles  which  are  Cu,  Ag,  Al2  O3 ,  CuO,  and  Ti02 .  It 
was  noticed  that  at  the  top  wall  and  for  Re  =  200,  an  enhancement  in 
the  Nusselt  number  was  registered  by  increasing  the  volume  frac¬ 
tion  of  the  nanoparticles.  Moreover,  there  was  an  enhancement 
in  Nusselt  number  at  the  top  wall  for  Re  >300  except  in  the  sec¬ 
ondary  recirculation  zone  where  insignificant  enhancement  was 


Table  2 

Summary  of  the  parameters  that  affect  the  flow  over  backward-facing  step. 
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registered.  It  was  found  that  the  high  Nusselt  number  inside  the 
recirculation  zone  is  mainly  depended  on  the  thermophysical  prop¬ 
erties  of  the  nanoparticles  and  it  is  independent  on  of  Reynolds 
number.  Both  the  Reynolds  number  and  thermophysical  proper¬ 
ties  of  the  nanoparticles  have  an  influence  on  the  value  of  Nusselt 
number  outside  the  recirculation  region.  However,  it  was  noticed 
that  outside  the  recirculation  zone,  nanoparticles  having  high  ther¬ 
mal  conductivity  (such  as  Ag  or  Cu)  have  more  enhancements  on 
the  value  of  the  Nusselt  number,  and  verse  visa  in  the  primary  and 
secondary  recirculation  zones. 

Very  recently,  the  same  authors  [170-172]  investigated  numer¬ 
ically  the  laminar  forced  and  mixed  convection  assisting  flow  over 
horizontal  and  vertical  backward-facing  steps  in  a  duct  using  differ¬ 
ent  nanofluids  (such  as  Au,  Ag,  AI2O3,  Cu,  CuO,  diamond,  SiC>2,  and 
Ti02).  It  was  found  that  a  primary  recirculation  region  developed 
after  the  sudden  expansion  and  it  starts  to  become  fully  developed 
downstream  of  the  reattachment  point.  The  reattachment  point  is 
found  to  move  downstream  far  from  the  step  as  Reynolds  number 
increases.  Moreover,  they  also  inferred  that  low  density  nanofluids 
have  higher  absolute  velocity  compared  to  high  density  nanoflu¬ 
ids.  It  is  found  that  a  recirculation  region  was  developed  straight 
behind  the  backward  facing  step  which  was  appeared  between  the 
edge  of  the  step  and  few  millimeters  before  the  corner  which  con¬ 
nects  the  step  and  the  downstream  wall.  In  a  few  millimeters  zone 
between  the  recirculation  region  and  the  downstream  wall;  a  U- 
turn  flow  was  developed  opposite  to  the  recirculation  flow  which 
is  mixed  with  the  unrecirculated  flow  and  travels  along  the  chan¬ 
nel.  Two  maximum  and  one  minimum  peaks  in  Nusselt  number 
were  observed  along  the  heated  downstream  wall.  It  is  inferred 
that  Au  nanofluid  has  the  highest  maximum  peak  of  Nusselt  num¬ 
ber,  while  diamond  nanofluid  has  the  highest  minimum  peak  in  the 
recirculation  region.  Nanofluids  with  higher  Prandtl  number  have 
higher  maximum  peak  of  Nusselt  numbers  after  the  separation  and 
recirculation  flow  is  vanished. 
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11.  Conclusions 

In  this  article,  a  comprehensive  review  of  previous  efforts  is 
presented  for  different  convective  flow  regimes  and  heat  transfer 
through  a  duct  having  backward  facing  step.  The  effects  of  sev¬ 
eral  parameters  in  geometry,  boundary  conditions,  and  types  of 
fluids  were  extensively  introduced  and  investigated.  In  backward¬ 
facing  step,  the  Nusselt  number  increases  as  the  parameter  such 
as  step  height,  velocity,  Prandtl  number,  and  aspect  ratio  increase, 
and  a  vice  versa  for  expansion  ratio.  It  was  reported  that  the 
reattachment  length  increases  as  step  height,  velocity,  Prandtl 
number,  aspect  ratio,  and  expansion  ratio  less  than  2.25  increase, 
and  decreases  with  increasing  expansion  ratio  greater  than  2.25. 
These  features  are  summarized  in  Table  2.  In  addition,  an  extensive 
review  for  preparation,  parameters,  mechanisms,  characteristics, 
convective  heat  transfer  enhancement,  applications,  safety,  and 
environmental  impacts  for  nanofluids  was  also  reported.  It  is  found 
that  nanofluids  will  solve  many  technological  problems  and  space 
limitation,  and  will  be  used  commercially  such  as  other  ordinary 
coolants  in  the  future.  It  was  observed  from  the  literature  review 
that  many  work  have  been  done  using  conventional  fluids  to  study 
the  heat  transfer  and  fluid  flow  characteristics  in  BFS.  Thus,  more 
work  using  nanofluids  is  needed  to  investigate  the  heat  transfer 
enhancement  over  backward  facing  step. 
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